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It  was  shown  that  the  stationary  phase  approximation  of  the  full  wave  solu¬ 
tion  is  precisely  equal  to  the  physical  optics  solution  for  scattering  by  rough 
surfaces.  Thus,  if  the  principal  contributions  to  the  scattered  fields  come 
from  the  regions  around  specular  points  of  the  rough  surface,  at  high  frequen¬ 
cies,  the  physical  optics  approximations  are  valid.  On  the  other  hand,  if  the 
principal  contributions  to  the  scattered  fields  do  not  come  from  specular  points 
of  the  rough  surface,  the  physical  optics  solutions  are  not  valid  even  if  the 
surface  meets  the  radii  of  curvature  criteria  associated  with  the  Kirchhoff  ap¬ 
proximations  for  the  surface  fields.  Similarly  it  was  shown  that  if  the  scale 
of  the  surface  roughness  and  the  slopes  are  6mall,  the  full  wave  solutions  re¬ 
duce  precisely  to  the  perturbation  solution  for  the  scattered  fields.  Thus  the 
full  wave  solutions  which  account  ‘for  both  specular  point  scattering  as  well  as 
Bragg  scattering  in  a  uniform  self-consistent  manner,  resolve  the  apparent  dis¬ 
crepancies  between  the  physical  optics  and  perturbation  theories.  Since  at  near 
normal  incidence  the  backscattered  fields  are  primarily  due  to  specular  point 
scattering  while  the  backscattered  fields  near  grazing  incidence  are  primarily 
due  to  Bragg  scattering,  a  two-scale  model  "mostly  based  on  physical  considera¬ 
tions"  was  advanced.  These  composite  models  are  able  to  explain  features  in 
radar  cross  section  that  no  theory  can.  A  problem  which  is  directly  associated 
with  the  two-scile  model  of  the  composite  rough  surface  concerns  the  specifica¬ 
tion  of  the  wavenumber  k,  where  spectral  splitting  between  the  large  and  small 
scale  surfaces  is  assumed  to  occur.  Furthermore,  the  two-scale  model  of  the 
composite  surface  is  restricted  by  the  assumption  that  the  large  and  small  sur¬ 
face  heights  h£  and  hfi  are  statistically  independent.  In  order  to  examine  the 
earlier  results  based  on  the  two-scale  models  of  composite  surfaces,  the  full 
wave  approach  (which  accounts  for  both  specular  point  scattering  as  well  as 
Bragg  scattering  in  a  uniform  self-consistent  manner)  is  also  applied  to  the  two- 
scale  model.  Thus  in  Sections  2  and  3  of  this  report  the  composite  surface  is 
decomposed  into  a  large  and  small  scale  surface. 

In  Section  4  the  full  wave  approach  is  used  to  determine  the  scattering 
cross  sections  without  assuming  the  two-scale  model  of  the  surface.  Thus  in  thii 
case  the  problems  associated  with  the  specifications  of  k.  do  not  arise  and  it 
is  necessary  to  assume  that  the  large  and  small  scale  surfaces  are  statistically 
Independent.  Using  the  full  wave  approach  specular  point  scattering  and  Bragg 
scattering  are  not  artificially  separated  out  in  the  expressions  for  the  scatter¬ 
ing  cross  sections. 

It  is  shown  that  while  the  values  for  the  like  polarized  backscatter  cross 
sections  obtained  from  the  unified  and  two-scale  versions  of  the  full  wave  solu¬ 
tions  are  in  substantial  agreement  for  all  angles  of  incidence,  there  are  very 
significant  differences  between  the  values  of  the  cross  polarized  backscatter 
cross  sections  obtained  from  the  uniform  and  two-scale  models  of  the  rough  sur¬ 
face.  This  is  shown  to  be  primarily  because  the  physical  optics  solutions  for 
the  cross-polarized  backscatter  cross  sections  vanish.  However,  as  is  shown  I 
through  the  use  of  the  unified  full  wave  approach,  only  at  the  specular  points 
of  the  rough  surface  Is  there  no  depolarization  of  the  backscattered  fields. 

Thus  a  two-scale  model  based  on  a  perturbed-physlcal  optics  approach  should  not 
be  used  to  determine  the  cross  polarized  scattering  cross  sections. 

In  Section  5  a  summary  of  the  proposed  research  for  the  next  six  months  is 
presented  and  suggestions  for  future  investigations  are  also  listed. 
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1.0  INTRODUCTION 


In  the  first  Annual  Technical  Report  (March  1,  1981-February  28,  1982), 
the  principal  elements  of  the  full  wave  approach  and  its  relationship  to 
earlier  solutions  of  scattering  problems  were  summarized. 

It  was  shown  that  the  stationary  phase  approximation  of  the  full  wave 
solution  is  precisely  equal  to  the  physical  optics  solution  for  scattering 
by  rough  surfaces.  Thus  if  the  principal  contributions  to  the  scattered 

fields  come  from  the  regions  around  specular  points  of  the  rough  surface, 

2  2  2 
at  high  frequencies,  kQ<h  >  >>  1,  (kQ  -electromagnetic  wave  number,  <h  > 

-mean  square  height),  the  physical  optics  approximations  are  valid.  On 

the  other  hand  if  the  principal  contributions  to  the  scattered  fields  do 

not  come  from  specular  points  of  the  rough  surface,  the  physical  optics 

solutions  are  not  valid  even  if  the  surface  meets  the  radii  of  curvature 

criteria  associated  with  the  Kirchhoff  approximations  for  the  surface  fields. 

Similarly  it  was  shown  that  if  the  scale  of  the  surface  roughness  is 
2  2 

small  (kQ<h  >  «  1)  and  the  slopes  of  the  rough  surface  are  small,  the  full 
wave  solutions  reduce  precisely  to  the  perturbation  solution  for  the 
scattered  fields  (Rice,  1951). 

Thus  the  full  wave  solutions  which  account  for  both  specular  point 
scattering  as  well  as  Bragg  scattering  in  a  uniform  self-consistent  manner 
resolve  the  apparent  discrepancies  between  the  physical  optics  and 
perturbation  theories. 

In  an  attempt  to  account  for  the  fact  that  at  near  normal  incidence 
the  backscattered  fields  are  primarily  due  to  specular  point  scattering 
while  the  backscattered  fields  near  grazing  incidence  are  primarily  due  to 
Bragg  scattering,  a  two-scale  model  "mostly  based  on  physical  considerations" 
was  advanced  (Wright,  1968;  Valenzuela,  1968).  These  composite  models  are 
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able  to  explain  features  in  radar  cross  section  that  no  theory  can.  More 
recently  Brown  (1978)  used  a  combination  of  Burrows'  perturbation  theory 
and  physical  optics  to  show  that  the  scattering  cross  sections  for  rough 
surfaces  can  be  expressed  as  a  sum  of  two  cross  sections;  one  associated 
with  the  filtered  surface  consisting  of  the  large  scale  spectral  components 
of  the  surface  and  the  second  associated  with  the  surface  consisting  of 
the  small  scale  spectral  components.  The  results  of  Valenzuela  (based 
on  physical  considerations)  and  those  of  Brown  (based  on  a  perturbed 
physical  optics  approach)  are  in  agreement  provided  that  the  mean  square 
slope  of  the  rough  surface  is  small. 

On  applying  the  full  wave  approach  to  the  two-scale  model  of  the 
composite  surface  it  was  shown  that  the  scattering  cross  sections  can  be 
expressed  as  a  weighted  sum  of  two  cross  sections.  It  was  also  shown  that 
the  difference  between  Valenzuela's  solution  and  Brown's  solution  for  the 
scattering  cross  sections  is  primarily  due  to  the  fact  that  in  Valenzuela's 
work  the  small  scale  surface  height  is  measured  perpendicular  to  the 
large  scale  surface  and  that  the  correlation  distance  for  the  small  scale 
surface  depends  on  distances  measured  along  the  large  scale  surface. 

On  the  other  hand  in  Brown's  work  the  small  scale  surface  height  is 
measured  perpendicular  to  the  mean,  reference  surface  and  the  correlation 
distance  is  measured  in  the  mean  Reference  plane.  This  is  contrary  to 
Burrows'  perturbation  formulation  upon  which  Brown's  analysis  is  based. 

A  problem  which  is  directly  associated  with  the  two-scale  model  of 

the  composite  rough  surface  concerns  the  specification  of  the  wavenumber  k, 

d 

where  spectral  splitting  between  the  large  and  small  scale  surfaces  is 
assumed  to  occur  (Brown,  1978).  Using  Brown's  approach  k^  is  specified 


on  the  basis  of  the  features  of  the  small  scale  surface  h  :  namely 

s 

2  2 

$  *  4ko<hg>.  It  is  shown  that  Brown's  results  for  the  scattering  cross 
sections  critically  depend  on  the  value  chosen  for  B.  Since,  using  Brown's 
theory,  the  small  scale  surface  must  satisfy  the  perturbation  criteria, 

Brown  suggests  that  k^  should  be  chosen  subject  to  the  condition  6  =  0.1. 

On  the  other  hand  Hagfors  (1966)  and  Tyler  (1976)  suggest  that  for  near  normal 
incidence  the  specification  of  k^  should  be  based  on  the  characteristics 
of  the  large  scale  surface  since  the  large  scale  surface  must  satisfy  the 
radii  of  curvature  criteria  associated  with  the  Kirchhoff  approximations 
of  the  surface  fields. 

In  addition  to  the  above  problems  associated  with  the  specification 
of  k^,  it  should  be  noted  that  the  two-scale  model  of  the  composite 
surface  is  restricted  by  the  assumption  that  the  large  and  small  surface 

heights  h  and  h  are  statistically  independent  (Brown,  1978) .  In  order 

jL  s 

to  examine  the  earlier  results  based  on  the  two-scale  models  of  composite 
surfaces  the  full  wave  approach  (which  accounts  for  both  specular  point 
scattering  as  well  as  Bragg  scattering  in  a  uniform  self-consistent  manner) 
is  also  applied  to  the  two-scale  model.  Thus  in  Sections  2  and  3  of  this 
report  the  composite  surface  is  decomposed  into  a  large  and  small  scale 
surface.  However,  the  value  of  B  (which  directly  specifies  k,)  is  allowed 
to  vary  from  0.1  to  2.0  since  use  of  the  full  wave  approach  imposes  no 
restriction  on  the  specification  of  k^.  Thus  assuming  that  the  surface 
height  spectral  density  function  (or  its  Fourier  transform  the  surface 
height  autocorrelation  function)  is  known,  the  value  of  k^  (where  spectral 
splitting  is  assumed  to  occur)  is  determined  by  its  relationship  to 
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2  2 

8  *  4k  <h  >  (see  Section  3).  The  portion  of  the  spectral  density  function 
o  s 

W(k)  corresponding  to  the  filtered  large  scale  surface  h^  (k  <  k^)  is 
used  to  determine  the  mean  square  slope  of  the  large  scale  surface  h^.. 

This  in  turn  determines  the  probability  density  function  for  the  large 
scale  surface  slopes  since  the  random  surface  height  is  assumed  to  comprise 
a  superposition  of  a  sufficiently  large  number  of  zero  mean  independent 
component  heights  such  that  the  surface  height  and  all  of  its  derivatives 
are  Gaussian  (Brown  1978).  Thus  all  the  data  needed  to  compute  the 
scattering  cross  sections  (based  on  the  two-scale  model)  can  be  obtained 
from  the  surface  height  spectral  density  function  or  the  corresponding 
surface  height  autocorrelation  function  (see  Sections  2  and  3).  It  is 
shown  that  for  values  of  8  between  one  and  two  the  scattering  cross 
sections  based  on  the  full  wave  solutions  are  practically  insensitive  to 
variations  in  k^.  For  these  values  of  8  the  filtered  surface  h^  satisfies 
the  radii  of  curvature  criteria  as  well  as  the  criteria  for  deep  phase 
modulation  (see  Sections  2  and  3).  However,  the  perturbed-physical  optics 
approach  cannot  be  used  for  8  >  0.1  (Brown,  1978). 

In  Section  4  the  full  wave  approach  is  used  to  determine  the  scattering 
cross  sections  without  assuming  the  two-scale  model  of  the  surface.  Thus 
in  this  case  the  problems  associated  with  the  specifications  of  k^  do  not 
arise  and  it  is  not  necessary  to  assume  that  the  large  and  small  scale 
surfaces  are  statistically  independent  (Brown,  1978).  Thus  using  the  full 
wave  approach  specular  point  scattering  and  Bragg  scattering  are  not 
artificially  separated  out  in  the  expressions  for  the  scattering  cross 


sections. 
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It  is  shown  in  Section  4  that  while  the  values  for  the  like  polarized 
backscatter  cross  sections  (<a  >and  <o  >) obtained  from  the  unified  and  two 
scale  versions  of  the  full  wave  solutions  are  in  substantial  agreement  for 
all  angles  of  incidence,  there  are  vary  significant  differences  between  the 
values  of  the  cross  polarized  backscatter  cross  sections  obtained  from  the 
unified  and  two-scale  models  of  the  rough  surface.  This  is  shown  to  be 
primarily  because  the  physical  optics  solutions  for  the  cross-polarized  back¬ 
scatter  cross  sections  (associated  with  the  filtered  surface  h^)  vanish  (see 
Section  4).  However,  as  is  shown  through  the  use  of  the  unified  full  wave 
approach,  only  at  the  specular  points  of  the  rough  surface  is  there 
depolarization  of  the  backscattered  fields.  Thus  no  matter  how  one  chooses 
kj  (the  wavenumber  where  spectral  splitting  is  assumed  to  occur)  a  two 
scale  model  based  on  a  perturbed-physical  optics  approach  should  not  be 
used  to  determine  the  cross  polarized  scattering  cross  sections. 

For  the  convenience  of  the  reader  of  this  report,  all  the  principal 
analj.i.cal  results  leading  to  full  wave  expressions  for  the  normalized 
cross  sections  are  presented  in  Section  (1 . 1) .  Thus,  starting  with  Maxwell's 
equations  for  the  transverse  components  of  the  electromagnetic  fields,  the 
boundary  conditions,  the  expressions  for  the  complete  field  expansions  and 
the  associated  orthogonal  relations,  the  generalized  telegraphists  equa¬ 
tions  are  derived.  It  is  shown  how  these  coupled  equations  are  solved 
using  a  second  order  iterative  approach  and  the  judicious  use  of  a  local 
coordinate  system  that  conforms  with  features  of  the  rough  surface.  Expli¬ 
cit  expressions  for  the  scattered  radiation  fields  and  the  like  and  cross 
polarized  scattering  cross  sections  are  also  presented  in  this  section. 


Thus  the  reader  finds  in  this  section  a  convenient  review  of  all  the  back¬ 
ground  material  that  leads  up  to  the  general  full  wave  expressions  for  the 
scattering  cross  sections. 

In  Section  (1 .2)  the  technical  reports  and  publications  resulting  from 
this  Air  Force  Contract  are  listed.  In  Section  5  a  summary  of  the  pro¬ 
posed  research  for  the  next  six  months  is  presented  and  suggestions  for 
future  investigations  are  also  listed. 
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1 . 1  Summary  of  Research 

In  this  section  the  principal  steps  leading  to  the  derivation  of 
the  full  wave  solutions  for  the  scattering  cross  sections  for  random 
rough  surfaces  are  summarized.  The  complex  electromagnetic  parameters 
of  the  media  above  and  below  the  irregular  boundary  y  =  h(x)  are  (see 


Fig.  1.1 
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£q»  y  >  h(x) 


e^,  y  <  h(x) 


(1.1) 


in  which  e^,  e^,  and  the  conductivity  a  are  real  and  an  exp(iiot)  time 
dependence  is  assumed. 

In  this  work  it  is  assumed  that  the  complex  permittivity  e  and 
permeability  y  are  not  functions  of  x.  Thus  wave  scattering  is  only  due 
to  the  irregular  boundary  between  the  two  half  spaces  y  >  h(x)  (medium  0) 
and  y  <  h(x)  (medium  1).  Since  it  is  convenient  to  represent  electric 
current  loops  by  magnetic  dipoles,  both  electric  and  magnetic  sources  J 
and  M  are  assumed  to  be  present.  Maxwell's  equations  for  the  transverse 
electric  and  magnetic  fields,  and  respectively,  are  (Bahar  1973a, b) 
3E 

--£•  *  \)  -  1ST  Vt'<Rt  *  ®x> 


+  M_  x  a  +  - —  v  j 
T  x  iioe  T  x 


(1.2) 


-g  -  *  v  - 1^7  vT-<s  x  st> 


+  a  x  +  T—  V  M  , 
x  T  iojy  T  x 


(1.3) 


in  which  the  operator  VT  is  given  by 

„  -  3  -  3 

V_  *  a  - —  +  a  - — 

T  y  3y  z  3z 


(1.4) 


>>>»;.' 


■w:>: 


.*  V  V  ^  ' 


Planes  of  incidence  and  scatter  with  respect  to  the  reference  coordinate 
system.  Mean  (reference)  plane  for  rough  surface  is  y  =  0. 


9 


! 


n 


m 

1 


and  the  transverse  vectors  are 


(1.5) 


A_  =  a  A  +  a  A  ,  A  =  E.H.J  or  M 
T  y  y  z  z  * 

The  following  field  transform  pairs  provide  the  basis  for  the  complete 

expansion  of  the  transverse  electric  and  magnetic  fields  into  vertically 

and  horizontally  polarized  waves: 


ET(x,y,z)  *  1  [E  (x,v,w)eT  +  E  (x,v,w)eT]dw, 


(1.6) 


where 


EP(x,v,w)  -  /"  ET(x,y,z)« (h^  -a^dydz,  P=V  or  H, 


HT(x,y,z)  =  £  /“  [HV(x,v,w)h^,  +  HH(x,v,w)h^]dw, 


rV  .  „H, 


.  rHi 


(1.7) 

(1.8) 


where 


HP(x,v,w)  *  /"  HT(x,y,z) • (a^  x  eP)dydz,  P=V  or  H, 


(1.9) 


The  basis  functions  for  the  vertically  polarized  waves  are 
iJJ  -  ZV(a  i|>V(v ,y)  -  - - H(w,z) 


u  +  w 


(1.10) 


and 


h^  -  az\J;V(v,y)i))(w,z) 


(1.11) 

and  the  complementary  basis  functions  for  the  vertically  polarized  waves 


are 


-  zV(l/(v,y)  +  2 


(1.12) 


and 


hy  =  azN>V(v,y)$C(w,z)  .  (1.13) 

For  the  horizontally  polarized  waves,  the  basis  functions  and  the  complementary 
basis  functions  are  respectively 

(1.14) 


e^  **  az^H(v,y)<*i(w,z)  , 


y  *  ’ ' >  'V'  «.** •  ' 


•  .V,'- V*  ,'«V« 


Rio'v  "  Ricrccr 


0  P 

C  +  C1K 


k0C0  ’  V1  *  klCl 


ni/n0  » 

n0/ni  ’ 


P=H  . 


(1.23) 


The  Intrinsic  impedance  for  medium  m  is  q  =  (p  /e  )  .  Furthermore 

m  mm 

in  (1.19)  and  (1.20) 


1  +  4and  Ipl-  1  +  njj 


(1.24) 


are  the  transmission  coefficients  and  th  (v,h)  is  the  value  of  the  surface 

s 

wave  basis  function  at  the  boundary  y  =  h.  For  a  homogeneous  model  of 
conducting  ground  the  surface  wave  exists  for  the  vertically  polarized 
waves  only.  The  superscripts  or  subscripts  V  and  H  are  used  to  denote 
vertically  or  horizontally  polarized  waves.  The  symbol  Zv  denotes  the 
summation  (integration)  over  the  entire  wavenumber  spectrum  v  consisting 


of  the  radiation  term  (0  <  Vq  <  <*>) ,  the  lateral  wave  term  (0  <  v^  <  °°) , 


and  the  guided  surface  wave  term  v^.  For  i  *  0  or  1 
v±  *  (k2  -  u2  -  w2)1*,  Im(vi)  <  0 


(1.25) 


and  for  the  surface  wave  term  v^g  is  a  solution  of  the  modal  equation 


-1/R^0  -  1/B^  -  0  P-V  or  H 
The  wavenumber  for  medium  i  =  0,  1  is 


(1.26) 


ki  =  “^iei^ 


(1.27) 


The  basis  functions  satisfy  the  biorthogonal  relationships  for  P  and  Q  equal 
to  V  or  H 


eT  • (hQ  x  ax) 'dydz 


hjji  • (ax  x  e£) 'dydz 


-  6p  qA (v-v ' ) A (w-w ' ) 


(1.28) 


in  which,  for  the  primed  quantities  the  variables  are  u',v',  and  w'  and 
A (v-v')  =6 


6  (v-v'),  vVv 


6 

V,V  ,  v'-V 


(1.29) 


The  subscripts  q,r  are  equal  to  0,  1  or  s  for  the  rad iation, lateral  wave  and 
surface  wave  terms,  respectively,  and  6(a-8)  and  6  are  the  Dirac  delta 

Qt  9  P 

function  and  Kronecker  delta,  respectively.  The  completeness  relationship  is 


6(y-y0)  *  Z  I?  ^P(v,y)<|<P(v,y0) 


!p  n l  ^(v«y^5(v*yo)dvo 


+  IP  NP  4'P(v,y)4'P(v,yQ)dv1 


+  IP  NP  4'P(v,y)^P(v,yQ) 


(1.30) 


The  scalar  basis  functions  for  the  radiation  fields,  and  the  lateral  and  surface 
P  P  P 

waves  are  ipn,  ip  ,  and  ,  respectively.  The  corresponding  normalization 

U  1  8 

coefficients  are 


NP(v)  *  1,  NP(v)  “  R^j/2wIP(v),  j-0,1. 


(1.31) 


in  which  I,  are  the  transverse  wave  impedances  or  admittances.  For  j*0,l 
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I 


\ 


C 

i 

'K 


t 


> 


S 


s 


V  2  .  2 

zY(v)  - 

j  UUEj 


'j 


"2  -  ■  -  Zq(v) 


2  2  ,, 
u— -  Z?<v> 

UWE,  1 


*>) 


2,2 

u  +  w 


uoip. 


2  2  u 

u  .±^..  «  YH(v) 
uam  X(P  ' 


0 


2  2 

u  +  w  „H,  v 

— - - *  Y,  (v) 

uuy,  1 


(1.32) 


Thus  in  the  above  expressions  for  ijj  (which  are  independent  of  y) ,  the 


parameters  for  medium  0  are  used  for  the  radiation  term,  and  the  parameters 

.P. 


for  medium  1  are  used  for  the  lateral  wave  term.  On  the  other  hand  I  (v,y) 
which  is  a  function  of  y  is  given  by 


I  (v,y)  -  \ 


ZV(v,y)  - 


zo  .  y  >  h 


z1  ,  y  <  h 


YH(v,y) 


,H 


Y"  ,  y  >  h 


.H 


Y“  ,  y  <  h 


(1.33) 


The  irregular  boundary  and  its  gradient  are  assumed  here  to  be  continuous 
functions  of  x  only.  Thus  the  exact  boundary  conditions  at  y  *  h(x)  can  be 
expressed  exclusively  in  terms  of  the  transverse  field  components 


ISF  V  (5t  x  ;x> 


-E  £ 
y  dx 


■•HI*- 


(1.34) 


~h+ 


ln/~  r;  v  u  dh  “  Llh+  n 
iuu  T  x  T  y  dx  _  *  [_  zj  _ 


(1.35) 


|h  ”  h 

The  complete  field  expansions  are  substituted  into  Maxwell's  equations 
for  the  transverse  field  components  and  use  is  made  of  the  orthogonality 


...  .v-v^vjs . 


\-\  «\  .  -*.v.v  . . 


relationship.  Green's  theorem,  and  the  exact  boundary  conditions,  to  obtain 

P  P 

the  differential  equations  for  the  field  transforms  E  and  H  .  These  may 

P  P 

be  expressed  in  terms  of  the  forward  and  backward  wave  amplitudes  a  and  b 


respectively,  as  follows: 


„P  P  ^  ,  P  ,  ^P  P  -  ,  P 
H  ■  a  +  b  and  E  ■  a  +  b  ,  P 


V,  (upper  sign 


H,  lower  sign 


(1.36) 


Thus  Maxwell's  equations  are  converted  into  the  following  generalized 


telegraphists  equations  for  P«V  or  H  (Bahar,  1973a, b) 

.  Js!  .  luaf  .  j  £,/(sm  aQ  +  „Q>du.  .  Ar  , 


(1.37a) 


^  -  i-b-  -  5  ;./<«“  .« ♦  «£  b<W'  +  Bp  . 


(1.37b) 


The  transmission  scattering  coefficients  for  the  radiation  fields  are  given 
by  (Bahar,  1973b) 


S>’ 


i  R><v)  p  p 

,v')  1  0 - c;(v’,v)<(v,v' 


where 


2  p 


/expCi (w-w')z]dz  s  I(w,w')  -*■  6(w-w') 
-£ 


)  1 (w , w ' )  P“V ,H,  a H  (1.38a) 


for  £  -*>  00  and 

(v) 

SHV^v,v'^  “  “  SVH^v'»v^  ~\T -  m  ~  I  Ch^v‘v'^  I(w*w')  ( 

N0(v’) 

The  reflection  scattering  coefficients  for  the  radiation  fields  are 


(1.39) 


Nn<v) 


(v,v')  -y  - - C^(v',v)  +  C*(v,v')  I(w,w') 


pp .  2  pv  ’  '  p' 


N?(V)  .  y 

SHV(V,V,)  "  sVH(v',v)  v  ,  "  2  CH(v»v,)  1<w»w') 

N(PV  ) 


(1.40) 


(1.41) 


The  coupling  coefficients  Cp  are  given  by 


-»  *o(v-y)  af * 0(v’*y)) 

Jh” 


(1.42 


On  employing  the  properties  of  the  scalar  basis  function  it  can 


be  shown  that  (1.42)  reduces  to 


Cp(v',v) 


IP(v,y)N?(v') 


,2  2 

v  -  v 


v'2  ^(v.yty^v'  ,y) 


+  «'o(v-y)  a5'*o(''’iy>)  +  • 


(1.43 


Similarly,  it  can  be  shown  that 

T  iw»J(v')  [  y 

ch<v’-v>  ■  -w5 —  •  *o<v-y>  *o(v’-y) 


+  7  sf  *0<v-y)  'ta5?  *o(v'  •y) 

+  ihtuf+wU - 3  ,y(V(y)  i(.?(v',y) 


Thus  on  using  the  properties  of  the  scalar  basis  function  i|>q, 

.  „  luNo(y,)  dh 

CH(V-V>~  — 

-  ^y  *J(v.5r)l,‘o<v' ,y)  "  *0(',,y) 


.  -1*" 


9y  r0 


*„(v',y) 


W^n ( v ' )  Ah  V  H 

-V-i  *j(v.h)*;(v,,h) 

Jy.  fi .  3»]  - .  fx  -  fa  fl 

L1  l  “lJ  M  C1  JJ 


(1.44 


(1.45 


Since  excitations  of  vertically  and  horizontally  polarized  waves 

(with  respect  to  the  reference  (x,z)  plane)  are  considered,  both  vertical 

electric  and  magnetic  dipoles  located  at  r  ■  r*  ■  x^a  +  y^a  +  z^a  ■  -r*n* 

x  y  z 

are  treated  in  this  analysis.  Thus  for  vertical  electric  and  magnetic  dipoles  of 

P  P 

intensity  J  and  M  respectively  the  source  terms  A  and  B  appearing  in  (1.37) 
are  given  by 

A  (v)  ■  -  B  (v) 

“  -JZQ(v)NQ(v)i|/Q(v,yi) 

•  exp(iwz^)6  (x-x^)/4ir  ,  (1.46) 

AH(v)  «*  BH(v) 

-  MYQ(v)NQ(v)^Q(v,y1) 

•  exp(iwz*)6  (x-x*)/4tt  (1*47) 

P  P 

The  first-order  iterative  solutions  for  the  wave  amplitudes  a  and  b 
are  obtained  by  neglecting  the  transmission  and  reflection  scattering 
coefficients  in  (1.37).  These  first-order  solutions  are  substituted  on 
the  right  side  of  (1.37),  and  the  resulting  equations  are  solved  to  obtain 
the  second-order  iterative  solution  for  the  wave  amplitudes.  These 
second-order  iterative  solutions  are  used  in  the  complete  expansions 
for  the  electromagnetic  fields  to  obtain  the  desired  iterative  solutions 
for  the  scattered  radiation  fields  through  the  use  of  the  steepest 
descent  method.  Thus  the  first-order  solutions  to  (1.37)  are  the  unperturbed 
vertically  and  horizontally  polarized  fields  excited  by  the  vertical  electric 
and  magnetic  dipoles  respectively.  The  second-order  iterative  solutions 
which  account  for  depolarization  and  scattering  in  arbitrary  directions 
are  suitable  when  multiple  scattering  can  be  ignored.  Thus  for  x  >  x* 
the  unperturbed  wave  amplitude  excited  by  a  vertical  electric  dipole  of 
intensity  J  (amp-meter)  is 


a  (x,v,w)  -  H  (x,v,w) 


exp(-iux)  /^exp(iux,)AV(x')dx' 
exp[~iu(x-x1) ]Aq  . 


where 


.V  -J  r.,V  .V  - 
A0  "  4 it  ^0  '  r*r* 


(1.48) 

(1.49) 


and 


i>(w,z)  -  exp(+iwz)  (1.50) 

Substitute  (1.48)  into  (1.8)  to  obtain  the  unperturbed  magnetic  field 


HT(x,y,z).  Thus 


Hz(x,y,z) 


V  v  - 

H  (x,v,wH0(v,yH(w,z)dv0dw. 


(1.51) 


0 


Noting  that 


(1.52) 


it  can  be  shown  that  for  y  >  h  and  y  >  h,  (1.51)  can  be  written  as 


Hz(x,y,z)  -  Rvo*o  exP[~  i(ux  +  wz)l 


8tt 


exp[i(ux1'  +  wz*  -  v^y'^jjdVgdw  . 


(1.53) 


When  the  vertical  electric  dipole  is  far  from  the  rough  surface  (kr  >>  1) 
the  steepest  descent  method  is  used  to  evaluate  (1.53)  (Brekhovskikh,  1960) 


u  -  kosin0o  cos$ i  vQ  -  kQcose0, 


w  ■  kQsin0Q  sin$ 


(1.54a) 


i  i  .  Qi  ,  i  i  i  1 
x  -  -r  sin0Q  cos<}>  ,  y  -  r  cos0Q, 


and  z^  ■  -r^sin0Q  sin$^ 


(1.54b) 


and 


x  -  r  sin0^  cos$^ ,  y  -  r  cos0^  ,  and 


z  ■  r  sin©!  sin<^  . 


(1.54c) 


It  can  be  shown  that  the  unperturbed  magnetic  field  near  the  boundary  y  *  h, 
is  given  by 

H  (x,y,z)  =  exp(-ik  r1) 

2  Airr1  0 

•  [uR^J  exp  -  i(ux  +  wz)]  ±  (1.55) 

in  which  0Q  and  <p  are  replaced  by  0q  and  (ji1  in  the  expressions  for  u,  v, 
and  w,  and  is  the  direction  of  the  incident  waves  from  the  vertical 

dipole.  Hence  the  azimuthal  component  of  the  incident  unperturbed  magnetic 
field  at  the  origin  is 

H*  «*  -  H^/cos$  *  ikgJ  sin0g  exp(-ikgri)/4irr^'.  (1.56) 

To  obtain  the  radiation  fie Id » transform  use  is  made  of  (1.9).  On 
integrating  with  respect  to  y  and  z  and  using  the  biorthogonal  relationship 
(1.28)  we  obtain 

HV(x,v,w)  -  H^J0(vi)exp(-iuix) 

<5(w-w')6(v-v*-) 

=  aj  exp(-iuix)6(w-wi)6(v-vi)  -  aV  (x,v,w)  (1.57) 

where  2 l  is  the  dimension  of  the  surface  in  the  z-direction  (-&  <  z  <  i) . 

Thus  for  a  vertical  electric  dipole  at  a  large  distance  from  the  rough  surface 

the  unperturbed  field  transform  is  represented  by  a  Dirac  delta  function. 

To  obtain  the  second-order  iterative  solution  that  accounts  for  wave 

scattering, (1.57)  is  substituted  on  the  right-hand  side  of  the  coupled 

V 

differential  equations  for  the  wave  amplitude  a  (x,v,w).  The  resulting  differ¬ 
ential  equation  for  the  vertically  polarized  forward  scattered  wave  amplitude 


is  integrated  to  yield 


a  (x,v,w)  -  -exp(-iux)  j  Sw(v,v  ) 


•  a*  exp  i(u-u^)x,dx' 


(1.58) 


A  similar  expression  can  be  derived  for  the  back  scattered  wave  amplitude. 

For  an  irregular  boundary  of  dimension  2L  in  the  x-direction,  y  =  h(x) 

(-L  <  x  <  L) ,  the  forward  scattered  wave  amplitude  can  be  expressed  as  follows 

aV(x,v,v)  -  -Mte.G1)  ggf-ffi 

u  -u 


Li 

’ 

^  exp  i[(u-ui)x  +  (vQ+VQ)h]dx  .  (1.59) 


Integrate  (1.59)  by  parts  to  obtain  the  total  (specular  and  non-specular) 
scattered  wave  amplitude  (assuming  that  k^L  >>  1), 

HV(x,v,w)  -  aV(x,v,w)  »  M(0,8i)I(fl,fli,h,L)2L 


•  exp  (-iux)/[(u+ui)  (Vq+Vq)  ] 


(1.60) 


where  I(£J,£]  ,h,L)  is  given  by 


Li 

I(fi  ,0±,h,L)  -  exp{ik0[(cj  +  CQ)h 


+  (Sq  cob<J>  -  Sq  cos<()i)x]}dx. 


(1.61) 


The  forward  scattered  radiation  field  is  obtained  by  substituting  (1.60)  into 
(1.51)  using  the  relationships  (1.52) 


Hz(x,y,z) 


HV(x,v,w) 


exp-i(Vgy  +  wz)dVp  dw 
00  00 

_  [  M(n.n1)i(n.oi.i 

'  _  _  (u+u1) (v.+vj) 


•  exp-i(ux  +  v_y  +  wz)dv_  dw 


n  .62^ 


Using  (1.54)  to  express  u,Vq  and  w  in  terms  of  0q  and  <J>,  and  x,y,  and  z  in 
terms  of  0q  and  the  steepest  descent  method  can  be  used  to  evaluate  (1.62) 
for  kgr  >>  1.  Thus  it  can  be  shown  that 


Hz(x,y,z) 

^  2iri  uf  exp(-jkr)M(flf  .fl*,h,L)2L 

/  f  .  *w  f  .  U 
r(u  +  u  ) (vQ  +  vQ) 


(1.63) 


The  solutions  for  the  other  scattered  fields  are  obtained  in  a  similar  manner. 

Thus  the  iterative  solutions  for  the  vertically  and  horizontally  polarized 

Pf  Pf  Pf 

scattered  radiation  fields  G  -  E  *  HqH  can  be  expressed  as  follows 


,Vf 


,Hf 


G0  I(5£,5if?a#A)c3 


fFW 


,HV 


r.VH 


.HH 


Lvil 

G 

_Hi 

G 

(1.64) 


Vf  Hf  Vf  Hf 

in  which  E  ,  E  and  H  ,  H  are  vertically  and  horizontally  polarized 


components  of  the  scattered  electric  and  magnetic  fields,  respectively,  and 

i 


-ikQ4£L 

- —  exp[-ikQr  ]  ,  Cq  -  cosB^ 


2irr 


(1.65) 


The  projection  of  the  rough  surface  in  the  (x,z)  plane  is  ■  4L£,  and  the 
observation  point  is  at  rf  -  rfn*\  The  expression  for  the  integral  I  in  (1.64) 


is 


£  L 

T,-f  -i  -  . .  1 

I(n  ,n  »rs,A)  -  ~ 

y  -£  -l 


•exp[ikg(n^-n*) *rs]dx  dz 


(1.66) 


in  which  n*  and  n^  are  unit  vectors  in  the  directions  of  the  incident  and  the 

scattered  wave  normals  in  medium  y  >  h,  and  r  is  the  distance  from  the  origin 

s 

to  points  on  the  rough  surface. 

For  two  dimensionally  rough  surfaces 

r  “  +  h(x,z)a  +  za  -  r  -  (y-h)I  «  r  -  fi 


(1.67a) 


^  -V 


V 


where 


f(x,y,z)  *  y  -  h(x,z)  =  0 


(1.67b) 


is  the  equation  for  the  rough  surface,  and  the  unit  vector  normal  to  the 
rough  surface  is 

n  -  Vf / | Vf | 

=  siity  cosSa  +  cosya  +  siny  sin6a 
x  y  z 

=  (-h  a  +  a  -ha  )/(h2  +  1  +  h2)^  (1. 

x  x  y  Z  Z  2  ' 

where 


(1.68a) 


9h  .  ,  9h  . 

—  =  h  and  —  =  h 
9x  x  9z  z 


(1.68b) 


In  order  to  extend  the  solution  (1.64)  to  problems  of  scattering  by  two- 
dimensionally  rough  surfaces  h(x,z)  with  arbitrary  slopes,  while  preserving 
the  relatively  simple  forms  of  these  iterative  solutions,  the  rough  surface 
is  regarded  as  a  continuum  of  elementary  surfaces  of  varying  slope  and 
height  rather  than  a  continuum  of  horizontal  elementary  surfaces  of  varying 
height.  The  contribution  to  the  total  scattered  field  from  an  inclined 
elementary  surface  is  obtained  from  (1.64)  after  making  the  appropriate 
coordinate  transformations.  In  this  way  the  restrictions  on  the  maximum  slope 
of  the  rough  surface  are  removed  while  preserving  all  the  advantages  of 
using  the  full-wave  approach. 

The  contribution  to  the  total  scattered  field  from  an  elementary  horizontal 


surfaces  at  r  -  rg  of  width  dx  dz  can  be  expressed  in  matrix  form  as 

dGf  -  Gq  dS  G1 


(1.69) 


in  which  the  constant  Gq  is  defined  in  (1.65)^G1  is  the  incident  field  at 
the  origin,  and 
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dS(x,z)  =  CqF  exp[ikg(n^-  n*)-rg]dx  dz/A 


=  CJF  dl(x,z)  (1.70) 

The  expression  corresponding  to  dS(x,z)  for  an  elementary  surface 

dA  =  ndxdz/(n*a^)  *  ndxdz/cosy  (1.71) 

is  given  by  dC.  In  matrix  notation  it  is  expressed  as  follows: 


* 

4 


cJVf^.H1)!1 


■exp[ik0(nf-ni) Tg]dA-n/Ay 


(1.72) 


In  (1.70)  the  angles  of  incidence  and  scatter  with  respect  to  the  reference 
plane  (x,z)  are  replaced  by  the  local  angles  of  incidence  and  scatter  with 
respect  to  the  local  tangent  plane.  The  the  elements  of  the  matrix  F 


(1.64)  are  (Bahar,  1981a, b) 


cinFW  2cJncJn[(yi:CjncJncos(4.in-^in)-sJnsJn)(l-l/er)+(l-yr)cos(<^rn-1j»ln)] 


„in„fn 
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An.  i  u  r  /  t-y  in.  f  n  / 1 1  n  .  m .  .  m.  r  n.  *  .  .  . .  .  / ,  i  n  ,  m .  i 

0  C0  L(erCi  C!  cos(0  -<P  )-SQ  SQ  )  (l-l/wr)+(l-er)cos(<(>  -<j>  )] 


.fn  ,in.  „in„fn 


fn  .in. 


cJV®. 


(cJn-K3*n/nr) (CQn+cJn/nr) (cJn+Con) 


C -sin(»f%in)2C^nnr[(l-l/cr)C^n-(l-l/ur)C^n] 
0  (C0n+nrCin)  (C0I1+Cin/nr)  (C0n4€0n> 


cinpVH  sin(4.f%in)2cJncJnnr[(l-l/Mr)cJn-(l-l/er)cJn] 

0  (cJn-K:Jn/nr)(CQn+T1rcJn)(cJn-+cJn) 

in  which  the  dimensionless  quantities  nr»^r»cr»  anc*  are  the  refractive 
index,  relative  intrinsic  impedance,  relative  permittivity,  and  relative 


(1.73) 


permeability,  respectively 
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(1.74) 


The  cosines  and  sines  of  the  angles  of  Incidence  and  scatter  (with  respect 


to  the  local  coordinate  system)  0j!jn  and  e^11  in  free  space,  y  >  h(x,z),  are 


.  1.2) 

,in  „in 

-i  - 

0  "  CO860 

■  -n  •  n 

,fn  „fn 
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-f  - 
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(1.75a) 
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b0  1 

|n^  x  n| 

(1.75b) 

The  sines  of  the  corresponding  angles  in  medium  1,  y  <  h(x,z),  are  given  by 
Snell's  law 


_in  .  .in 
-  sin01 


S>r 


_fn  .  „fn 
S1  =  sin©1 


S0n/”r 


(1.76) 


Thus 


C^n  -  cos0*n  -  [l-CS*")2]*5 


„fn  _fn  r,  ,_fn,.2-ii5 

C±  =  cosej^  -  [i-(s1  )  } 


(1.77) 


The  cosine  and  sine  of  the  angle  between  the  planes  of  incidence  and  scatter 
in  the  local  coordinate  system  are  given  by 

cos(<J>^n-4>*n)  *  an  •  a”  (1.78a) 

nf  Hi 

and 


sin($fn-(J>in)  -  [ 


-n  -n  — | 


Hf  Hi 


(1.78b) 


incidence  and  scatter  and  local  coordinate  system  with  unit 


aHf  =  (n  Xav)/In  xavl 


(/xn)/|/  xnl 


(1.79) 


aHi  "  X*yi 


(/  xn)  / 1  /  xn  | 


(1.80) 


The  vertically  and  horizontally  polarized  components  of  the  incident  and 
scattered  electric  and  magnetic  fields  with  respect  to  the  local  plane  normal 
to  the  unit  vector  n,  are  denoted  by  the  subscript  n.  They  are  related  to 
the  components  with  respect  to  the  reference  plane  through  the  following 
transformations 


Si  EV1 


c1  EHi 


(1.81) 


,Hf  -  TfGfn 


-Sf  EVf 
\p  n 

_f  Hf 
C,  E 


(1.82) 


in  which  C*  and  S*  are  the  cosine  and  sine  of  the  angle  between  the  local 
fp  rl> 

plane  of  Incidence  and  reference  plane  of  incidence  normal  to  the  unit  vectors 

a^  and  a ,  respectively.  Thus  they  can  be  expressed  in  terms  of  the  scalar 

product  and  the  scalar  triple  product 

_i  . i  -n 

c  =  cos *  -  aHi  .  aHi 


sf  -  sin/  -  I^H/^n1] 


(1.83) 


Similarly,  and  are  the  cosine  and  sine  of  the  angle  between  the  local 
plane  of  scatter  and  the  reference  plane  of  scatter  normal  to  unit  vectors 


a„,  and  a„c,  respectively.  Thus 


The  full  wave  solutions  for  the  scattered  wave  amplitudes  can  therefore  be 
expressed  as  follows  In  matrix  form 


Gf  *  Gn  f  C«nTfFTi 


exp[ikn(n^-n*)  •  r  ]U(r  )dA  •  nG3 


=  G()C(iIf,5i)Gi 

In  which  the  shadow  function  U(r  )  is 

s 


U(t8) 


a. 85 


1  Illuminated  and  visible  region 
0  nonilluminated  or  nonvisible  region 


(1.86 


The  normalized  scattering  cross  section  per  unit  area  for  rough  surfaces  is 


tshimaru,  1978) 


PQ  .  .  f .  2 1  Pf  1 2.r  i _Qi  1 2, 
j  v  -  4*(r  )  I E  |  /(AjE*  |  ] 


(1.87 


in  which  is  the  projection  of  the  area  of  the  rough  surface  A  on  the 
reference  plane  normal  to  a^.  Thus  for  P,Q-V,H 


/Q.^|cPQ|2-A[DPQ(;)DPQ*,;,) 
y  y 

•  <F  -J')]0<?  Wr)  dxdy 

^  0  s  8  s  s  (S'S  h5’-;  ) 

po  y  y 

where  D  is  the  element  of  the  matrix 


(1.88 


D  -  cJ^FT1 

and  the  symbol  *  denotes  the  complex  conjugate. 


(1.89 


For  the  incoherent,  diffuse  field  the  scattering  cross  section  for  unit 


cross-sectional  area  is  given  by  (Ishimaru,  1978) 
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«jPQ>  -  4ir(rf)2<|Epf  -  <Epf  >  |  2>/Ay  |  EQ±  | 2 


(1.90) 


in  which  the  symbol  <  >  denotes  the  statistical  average. 

Thus  the  full  wave  expressions  for  the  normalized  scattering  cross  section 
per  unit  area  is 
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1.2  Interim  Technical  Reports  and  Publications  in  Scientific /Technical 


Journals 

Preprints  of  the  following  manuscript  "Scattering  Cross  Sections 
for  Composite  Models  of  Non-Gaussian  Rough  Surfaces  for  Which  Decorre¬ 
lation  Implies  Statistical  Independence"  were  submitted  to  the  contract 
monitor  for  publication  as  an  Interim  Technical  Report  (March  1,  1982- 
Nov.  30,  1982). 

Abstract  of  Interim  Report 


The  full  wave  approach  is  used  to  determine  the  scattering  cross 
sections  for  composite  models  of  non-Gaussian  rough  surfaces.  It  is 
assumed  in  this  work  that  the  rough  surface  heights  become  statistically 
independent  when  they  decorrelate,  thus  no  delta  function  type  specular 
term  appears  in  the  expression  for  the  scattered  fields.  The  broad 
family  of  non-Gaussian  surfaces  considered  range  in  the  limit  from  ex¬ 
ponential  to  Gaussian.  It  is  seen  that  for  small  angles  of  incidence, 
the  like  polarized  cross  sections  have  the  same  dependence  on  the 
special  form  of  the  surface  height  joint  probability  density,  but  for 
large  angles  the  scattering  cross  sections  for  the  horizontally  pola¬ 
rized  waves  are  much  more  sensitive  to  the  special  form  of  the  joint 
probability  density.  The  corresponding  results  for  the  depolarized 
backscatter  cross  section  are  also  presented.  The  shadow  functions 
are  shown  to  be  rather  insensitive  to  the  special  form  of  the  joint 
probability  density. 

Preprints  of  the  following  manuscripts  and  conference  papers 
were  submitted  to  the  Office  of  Public  Affairs  and  project  monitor 


for  review  and  approval  for  publication: 


"Propagation  of  Vertically  and  Horizontally  Polarized  Waves  Excited 
by  Distributions  of  Electric  and  Magnetic  Sources  in  Irregular 
Stratified  Spheroidal  Structures  of  Finite  Conductivity  Generalized 
Field  Transforms",  Canadian  Journal  of  Physics,  Vol.  61  No.  1, 
pp.  113-127,  1983. 

"Scattering  and  Depolarization  of  Electromagnetic  Waves  in  Irregular 
Stratified  Spheroidal  Structures  of  Finite  Conductivity — Full  Wave 
Analysis",  Canadian  Journal  of  Physics,  Vol.  61  No.  1,  pp.  128- 
139,  1983. 

"Scattering  Cross  Sections  for  Composite  Surfaces  That  Cannot  Be 
Treated  as  Perturbed-Physical  Optics  Problems."  Radio  Science,  in 
press. 

"Computations  of  Scattering  Cross  Sections  for  Composite  Surfaces 
and  the  Specification  of  the  Wavenumber  Where  Spectral  Splitting 
Occurs",  Submitted  for  review. 

"Scattering  and  Depolarization  by  Rough  Surfaces:  Full  Wave 
Approach",  Proceedings  of  the  SPIE  International  Technical  Symposium 
of  the  International  Society  of  Optical  Engineering,  Vol.  358  No.  28 
pp.  1-14,  August,  1982. 

6.  "Comparison  of  Backscatter  Cross  Sections  for  Composite  Rough  Sur¬ 
faces  with  Different  Mean  Square  Slopes",  International  Journal  of 
Remote  Sensing,  in  press. 

7.  "Shadowing  by  Non-Gaussian  Rough  Surfaces  for  Which  Decorrelation 
Implies  Statistical  Independence",  Radio  Science,  in  press. 

8.  Joint  International  IEEE/ APS  and  National  Radio  Science  Meeting  at 
the  University  of  New  Mexico,  May  24-28,  1982,  "Scattering  Cross 
Sections  for  Composite  Surfaces  with  Large  Mean  Square  Slopes — 

Full  Wave  Analysis." 

9.  The  SPIE  26th  Annual  International  Technical  Symposium  of  the 
International  Society  for  Optical  Engineering,  San  Diego,  Califor¬ 
nia,  Aug. 23-27,  1982.  Title  of  Invited  Paper,  "Scattering  and 
Depolarization  by  Rough  Surfaces." 

10.  International  IEEE/APS  Symposium  and  National  Radio  Science  Meeting 
at  the  University  of  Houston,  Texas,  May  23-26,  1983,  "Rough 
Surface  Scattering  that  Cannot  Be  Analyzed  Perturbed-Physical  Optics 
Approaches . " 

11.  International  IEEE/APS  Symposium  and  National  Radio  Science  Meeting 
at  the  University  of  Houston,  Texas,  May  23-26,  1983,  "Scattering 
Cross  Sections  for  Composite  Surfaces  and  the  Wavenumber  Where 
Spectral  Splitting  Occurs." 


2.0  SCATTERING  CROSS  SECTIONS  FOR  COMPOSITE  SURFACES  THAT  CANNOT  BE 


TREATED  AS  PERTURBED-PHYSICAL  OPTICS  PROBLEMS 

2.1  Background 

Perturbation  and  physical  optics  theories  have  traditionally  been 
used  to  derive  the  scattering  cross  sections  for  composite  surfaces  that 
can  be  regarded  as  small  scale  surface  perturbations  that  ride  on 
filtered,  large  scale  surfaces.  In  this  case  perturbation  theory  accounts 
for  Bragg  scattering,  while  physical  optics  theory  accounts  for  spec¬ 
ular  point  scattering.  However,  for  a  more  general  class  of  composite 
surfaces  that  cannot  be  decomposed  in  such  a  manner,  the  perturbed-physical 
optics  approach  cannot  be  used.  In  these  cases,  it  is  shown,  using  the 
full  wave  approach,  that  the  specular  scattering  associated  with  a  fil¬ 
tered  surface  (consisting  of  the  larger-scale  spectral  components),  is 
strongly  modified,  and  that  Bragg  scattering  and  specular  point  scatter¬ 
ing  begin  to  blend  with  each  other.  Since  the  full  wave  solution 
accounts  for  Bragg  scattering  as  well  as  specular  point  scattering  in 
a  self-consistent  manner,  it  is  not  necessary  to  filter  (decompose) 
the  composite  surface  to  evaluate  the  scattering  cross  sections  in  the 
general  case.  However,  filtering  the  composite  surface  enhances  one's 
physical  insight  as  to  the  validity  (or  lack  thereof)  of  the  perturbed- 
physical  optics  decomposition,  and  also  facilitates  the  numerical  eval¬ 
uation  of  the  scattering  cross  sections. 

2.2  Discussion 

In  order  to  account  for  Bragg  scattering  as  well  as  specular  point 
scattering  from  random  rough  surfaces,  composite  models  of  the  surface 
with  different  roughness  scales  have  been  considered  (Wright  1968, 


Barrick  1970,  Barrick  and  Peake  1968).  Thus,  for  instance,  for  back- 
scatter  more  than  30°  away  from  the  vertical,  it  is  shown  through  the 
use  of  perturbation  theory  (Rice  1951)  that  Bragg  scattering  is  domin¬ 
ant  and  the  scattering  cross  section,  which  is  polarization  dependent, 
is  proportional  to  the  surface  height  spectral  density  function.  On 
the  other  hand,  using  physical  optics  theory  (Beckmann  1968),  the  back 
scattering  cross  section  near  normal  incidence  is  shown  to  be  primarily 
due  to  specular  point  scattering  and  to  be  independent  of  polarization. 

Using  the  composite  models  of  Wright  (1968),  Semyonov  (1966)  and 
Valenzuela  (1968)  which  are  "mostly  based  on  physical  considerations," 
the  rough  surface  is  regarded  as  patches  of  slightly  rough  surfaces 
that  ride  over  the  large  waves.  Thus,  in  their  work  the  scattering 
cross  section  associated  with  the  surface  with  the  small  scale  roughness 
is  obtained  by  averaging  over  the  distribution  of  slopes  of  the  large 
scale  roughness,  or  by  averaging  over  the  tilt  angles  in  and  perpendic¬ 
ular  to  the  plane  of  incidence. 

More  recently  Brown  (1978,  1980)  applied  a  combination  of  Burrows' 
perturbation  theory  (1967)  and  physical  optics  theory  (Beckmann  1968) 
to  derive  the  backscattering  cross  section  from  a  perfectly  conducting 
two  scale  model  of  rough  surfaces.  The  first  term  in  his  solution  is 
the  specular  point  backscattering  cross  section  associated  with  the 
large  scale  surface  height  h^  and  the  second  term  is  the  Bragg  scatter¬ 
ing  cross  section  associated  with  the  small  scale  surface  height  h  . 

s 

Thus,  in  his  work  it  is  necessary  to  decompose  (i.e.,  spectrally  filter) 
the  composite  surface.  To  this  end.  Brown's  specification  of  the 


wavenumber  k^  (where  spectral  splitting  is  assumed  to  occur),  is  based 
entirely  upon  the  characteristics  of  the  small  scale  structure  (Brown 
1978).  However,  in  the  approaches  of  Hagfors  (1966)  and  Tyler  (1976) 
the  specifications  of  kj  (spectral  truncation)  is  assumed  to  be  based 
on  the  characteristics  of  the  large  scale  surface.  In  their  approaches, 
however,  they  ignore  the  effects  of  the  surface  comprising  the  small 
scale  spectral  components  (k  >  k^)  since  their  results  are  not  meant  to 
explain  scattering  far  removed  from  the  specular  direction.  In  addi¬ 
tion,  it  should  be  pointed  out  that  even  if  the  spectral  components  of 
the  large  scale  (filtered)  surface  satisfy  the  radii  of  curvature  cri¬ 
teria  (associated  with  the  Kirchhoff  approximation  of  the  surface 
fields),  they  may  not  necessarily  satisfy  the  condition  for  deep  phase 
modulation  implicit  in  the  evaluation  of  the  specular  point  result 
(Barrick  1970). 

In  order  to  apply  the  full  wave  approach  (Bahar  1981  a,b,  1982a, b) 

to  problems  of  scattering  by  rough  surfaces,  it  is  not  necessary  to 

decompose  (filter)  the  rough  surface  into  one  surface  h  with  a  small 

s 

roughness  scale,  and  another  with  large  radii  of  curvature.  However, 
when  such  a  decomposition  is  feasible  and  the  restrictions  on  both  the 
large  and  small  scale  roughness  are  satisfied  simultaneously,  the  full 
wave  solutions  for  the  scattering  cross  sections  (which  account  for 
both  Bragg  scattering  and  specular  point  scattering  in  a  self-consistent 
manner) ,  can  be  expressed  in  terms  of  a  weighted  sum  of  two  cross 
sections  in  agreement  with  perturbation  and  physical  optics  theories. 

A  comparison  between  the  full  wave  approach  and  the  approaches  of 
Brown  and  Valenzuela  has  also  been  made  recently  for  surfaces  with 


moderately  large  mean  square  slopes  (see  Section  2.3). 

The  major  objective  of  this  work  is  to  apply  the  full  wave  approach 
to  composite  surfaces  which  cannot  be  decomposed  into  a  two  scale  model 
without  violating  either  the  small  height  variance  criteria  essential 
to  the  application  of  perturbation  theory  or  the  large  radii  of  curva¬ 
ture  criteria  essential  to  the  application  of  physical  optics  theory. 
Thus,  in  Section  Q. .3)  the  principal  formulas  for  the  scattering  cross 
sections  are  developed  for  surfaces  that  do  not  satisfy  either  the  small 
variance  or  the  large  radii  of  curvature  restrictions.  In  Section (2.4) 
a  composite  surface  which  does  not  satisfy  the  physical  optics  (Kirch- 
hoff)  criteria  is  analyzed  in  detail.  To  decompose  this  surface  it  is 

assumed  that  the  wavenumber  k.  is  specified  in  accordance  with  the  cri- 

a 

teria  proposed  by  Tyler  (1976)  and  that  the  surface  that  rides  upon  the 
filtered  (Kirchhoff)  surface  may  not  satisfy  the  perturbation  criteria. 
It  is  also  assumed  that  if  the  wavenumber  kj  were  specified  in  accord¬ 
ance  with  the  criteria  based  on  the  variance  of  the  small  scale  surface 
(Brown  1978),  the  remaining  large  scale  surface  would  not  satisfy  the 
standard  physical  optics  criteria.  Thus,  a  perturbed-physical  optics 
approach  cannot  be  applied  to  the  problem.  It  is  shown  that  the  scat¬ 
tering  cross  sections  for  this  general  class  of  rough  surfaces  can  be 
expressed  as  a  weighted  sum  of  two  cross  sections.  The  contribution 
associated  with  the  filtered  (Kirchhoff)  surface  is  multiplied  by  a 
factor  which  is  significantly  smaller  than  unity  and  the  term  associated 
with  the  remaining  surface  that  rides  upon  the  filtered  surfaces  may 
contribute  significantly  even  in  the  near  specular  direction.  Thus  for 
this  general  model  it  is  illustrated  that  classic  Bragg  scattering  and 
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specular  point  scattering  begin  to  blend  with  each  other.  Since  the 
full  wave  approach  accounts  for  both  Bragg  and  specular  point  scattering 
it  is  not  necessary  to  decompose  (filter)  the  composite  surface  when 
the  full  wave  approach  is  applied  to  this  broad  class  of  scattering 
problems;  however  the  decomposition  may  provide  additional  physical 
insight  into  prior  approaches  and  facilitate  the  computation  of  the 
scattering  cross  sections. 

2.3  Formulation  of  the  Problem 

Using  the  full  wave  solutions  for  the  incoherent  scattered  radiation 
fields,  the  expression  for  the  normalized  scattering  cross  section  per 
unit  area  (Ishimaru  1978)  is  given  by 
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Second  order  iterative  solutions  for  the  scattered  wave  amplitudes 
are  used  in  the  derivation  of  (2.1)  (Bahar  1981a, b;  1982a, b).  This 
approximation  of  the  full  wave  solution  is  suitable  when  multiple  scat¬ 
tering  can  be  neglected.  In  (2.1)  Is  the  projection  of  the  rough 
surface  on  the  reference  plane  (y*0)  (see  Fig.2.1)and  kQ  Is  the  free 
space  wave  number  of  the  electromagnetic  wave.  The  vector  v  in  the 
cartesian  coordinate  system  (x,y,z)  is 

v  -  k^  -  k*  -  k  (n^  -n1)  «  v  a  +  v  a  +  v  a  (2.2) 

o  xx  yy  zz  '  # 

where  n*  and  n^  are  unit  vectors  in  the  direction  of  the  incident  and 
scattered  wave  normals  respectively.  The  coefficient  SP<^  is  (Bahar 
1981a) 
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Planes  of  Incidence  and  scatter 
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in  which  D  depends  explicitly  upon  the  polarization  of  the  incident 

wave,  (second  superscript  Q-V  -  vertical,  H  -  horizontal)  and  the 

polarization  of  the  scattered  wave  (first  superscript  P  “  V,H),  the 

directions  of  the  incident  and  scattered  wave  normals  n1  and  n^,  the 

complex  permittivity  and  permeability  of  the  medium  of  propagation  and 

the  unit  vectors,  n(h  ,h  ),  n^h'.h*)  normal  to  the  rough  surface  at 

X  Z  X  z 

r(x,h,z)  and  r ’ (x* ,h' ,z ') .  Thus  for  the  random  rough  surface 

f(x,y,z)  -  y-h(x,z)  -  0  (2.4a) 

Vf  *  7(y-h(x,z))  -  (-h  a  +  a  -  h  a  )  ■  nlVfl  (2.4b) 
y  xx  y  z  z  11 

in  which  the  components  of  the  gradient  of  h(x,z), 

hx  =  3h/3x,  hz  ■  3h/3z  (2.4c) 

are  random  variables.  The  shadow  function  U(r)  is  unity  when  the  sur¬ 
face  is  both  illuminated  by  the  source  and  visible  at  the  observation 
point  and  zero  otherwise  (Sancer  1969) .  When  the  surface  height  h  and 
slopes  (h  ,h  )  are  statistically  independent  (a  condition  that  holds  for 
Gaussian  surfaces  at  each  point) 


DP% 


dpqu 

exp(iv  h)>  -<  ^-~>x  (v  ) 
n*a  y  n*a  y 

y  y 


(2.5a) 


in  which  x(Vy) 

X  (v) 


y 


is  the  surface  height  characteristic  function 
■  <exp(iv  h)>  m  J'  exp(iv  h)p(h)dh 

y  y 


(2.5b) 


In  (2.5b)  p(h)  is  the  surface  height  probability  density  function.  The 
symbol  *  denotes  complex  conjugate  and  the  symbol  <>  denotes  the  statis¬ 
tical  average. t 

2  2 

+As  an  example,  for  a  Gaussian  surface,  x(Vy)  •  exp(-v^<h  >/2),  where 
2  2 

a  **  <h  >  is  the  mean-square  surface  height. 


If  at  high  frequencies  the  higher  order  derivatives  of  the  rough  sur¬ 
face  h  can  be  neglected  (Sancer  1969,  Bahar  1981b;  Bahar  1982a, b) 


<SPQ«p(lv 


(n*a  )  (n'*a  ) 

y  y 


exp{iVy[hx(x-x')+hz (z-z  * ) J } 


p(n,n* ,UU')dn  dn'd(UU’) 


(2.6) 


in  which  p(n,n’,UU')  is  the  joint  probability  density  function  of  the 
slopes  n(hx,hz),  n(hx,hz)  and  the  shadow  function  product  U(r)U(r')  and 


dn  =  dh  dh  ,  dn*  =  dh'dh' 
X  Z*  X  z 


(2.7) 


The  joint  probability  density  function  p(n,n,,UU')  can  be  expressed 
as  follows: 


p(n,nf,UU')  *>  p(n,n,)p(UU' |n,n’) 


(2.8) 


Assuming  that  h  and  h  are  independent  variables  (which  is  true  for 
X  z 

surfaces  with  isotropic  roughness) 


p(n.n')  -  p(h  ,h')p(h  ,h') 
X  X  z  z 


(2.9) 


For  a  rough  surface  height  with  a  Gaussian  distribution,  for  example. 


P«Vh'x) 


2tto  /(1-C 

X  X 


h  -  2C  h  h'  +  h' 

X  XXX  X 
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in  which  it  is  assumed  that 

2  2  2 
<h  >  “  0,  <h  >  *o  and  <h  h’>  ■  o  C 

X  X  X  XX  XX 


(2.10a) 


(2.  10b) 


For  the  homogeneous  isotropic  surface  assumed  in  this  work  the  normal¬ 
ized  autocorrelation  function  C  is  a  function  of  distance 

x 


r.  -  (x-x')a  +  (z-z*) a  -  x.a  +  za 
d  x  z  d  y  d  z 


(2.11) 


The  density  p(h  ,h’)  is  given  by  (2.10)  with  h  replaced  by  h  and 

Z  Z  X  z 


<h  > 

z 


0,  <h  >“a 
*  z  z 


a  ,  <h  h  >  ■  a  C  **  a  C 

X  Z  Z  Z  Z  XX 


The  two  point  conditional  density  function  p(UU'|n,n')  is  expressed  in 
a  form  similar  to  the  one  given  by  Sancer  (1969)  for  p(u|n).  Thus, 


pOJU'In.n')  =  P2(nt,n1  In.n^SOJU'-l)  +  [l-P2 (nf  .n1  j  ^.S'^fiOJU’),  (2.13) 
in  which  <5 (a) is  the  Dirac  delta  function  and  P2(n^,n*|n,n')  is  the 
probability  that  the  points  r  and  r'  on  the  rough  surface  will  be  both 
illuminated  by  the  source  (incident  wave  normal,  n^(0i,<j>'i') )  and  visible 
at  the  observation  point  (scatter  wave  normal,  n^(0^,<|>^))  given  the 
value  of  the  unit  vectors  normal  to  the  surface  n  and  n',  at  these 


points.  In  (2.13) 
P2 (nf ,ni|n,n')* 


P2(nf,ni|n)  ,  for  r  -+•  r'.C^  1 


P2 (n  ,n  jn)P2(n  ,n|n'),  for  |r-r'|-*-  °°,Cx  -*■  0  (2.14) 


and 


P2(nf,ni|n)  =  P2(nr,n1|ns)S(nI*n)S(-n1*n) 


-f  -i i - 


,-f  -v 


-i  -x 


(2.15) 


in  which  ?2  (n^,ni|n)  is  the  probability  that  a  point  on  the  rough  sur¬ 
face  is  both  illuminated  and  visible  given  the  value  of  the  slopes  at 

the  point  and  P„  (n^,n*|n  )  is  its  value  at  the  specular  points  (n-m  ) 
z  s  s 

(Smith  1967,  Sancer  1969).  The  arguments  of  the  unit  step  functions 
S(-ni*n)  and  S(n^*n)  vanish  at  points  of  the  rough  surface  where  the 
incident  and  scattered  waves  are  tangent  to  the  surface.  Thus  S(-ni*ns) 
=  1  and  S(n^»ns)  *  1  (Smith  1967,  Brown  1980).  In  view  of  (2.14)  and 


(2.15)  it  is  assumed  that 


r.  \  v 


2-(C2+C2)/2 

X  z 


,(nf ,nijn,n')  -  |p2(nf ,n1|ng) j 

•S(nf*n)S(nf*n,)S(-ni‘n)S(-ni*n') 


(2.16) 


Substitute  (2.16)  into  (2.6)  and  integrate  with  respect  to  UU'  to  get 

<SP^exp[iv  (h-h')]>*  f- — - il—L  exp[iv  (h  x.+h  z  )] 

y  J  (n.a)(n»*a.)  y  xd  zd 


P2(n^,n'*"|n,n,)p(n,n,)dn  dn' 


(2.17) 


Since  <S  ^exp[ivy (h-h' ) ]>  is  assumed  to  be  a  function  of  the  distance 
Ir^) ,  (2.12) ,  the  scattering  cross  section  (2.1)  reduces  to 

oo 

<aPQ>  =  —  j fsPQexp[ivy(h-h-)]>  -  D  Qp2(n  »n  ln)  x(Vy)  2| 


•exp[iv  x,  +  iv  z, ]  dx.dz, 
r  x  d  z  d  d  d 


(2.18) 


In  (2.18)  the  integrand  vanishes  when  |r^|  is  much  larger  than  the 

correlation  distances.  The  evaluation  of  the  scattering  cross  section 

(2.18)  simplifies  considerably  if  it  can  be  assumed  that  for  |rd|  less 

than  the  correlation  distances  &  =  Z  (C  (i  )  «  C  (2.  )  ■  1/e  and  e  is 

x  z  x  x  z  z 

the  Neperian  number)n(r)  =  n'(r').  In  this  case  p(n,n')  ■+■  p(n)6(n-n’) 


:PQ  JnPQ,r 


D  x(r)U(r) 


(2.19) 


Thus,  for  instance,  when  the  radii  of  curvature  of  the  rough  surface 
are  very  large  compared  to  the  electromagnetic  wavelength  X,  or  when  the 
slope  of  the  surface  is  small  n-  ay 

/)  PQ  -  2 

<SP^exp[iv  (h~h' )]>'-</  -- — exp[iv  (h-h')] 
y  ^  ri-  a  y 

y 

•p(hx,hz)p2(nf,ni|n)dhxdhz>  (2.20) 


Equations  (2.1)  and  (2.20)  can  be  applied  directly  to  composite 
surfaces  (Brown  1978,  1980)  given  by 


h  »  h.+  h 
l  s 


(2.21) 


where  h.  and  h  are  assumed  to  be  statistically  independent  random  func¬ 


tions 


k*  <hZ>  -  k2  a2  «  1,  — r—  -  ]h  ]«  l,\~~  «  h  i«  1, 

o  c  n  c  *  a-v  I  ov  ’  1  S7  * 
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(2.22a) 
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(2.22b) 

and  the  radii  of  curvature  of  the  surface  h^  is  much  larger  than  the 
wavelength  A.  In  this  case  the  full  wave  solution  for  the  scattering 
cross  section  accounts  for  both  specular  point  scattering  as  well  as 
Bragg  scattering,  and  (2.21)  reduces  to  (Bahar  1981b,  1982a, b) 


«,*«>  .  <„«>>  +  «,*!>, 

O  1 


(2.23) 


The  first  term  in  (2.23),  <oP^>  is  the  specular  scatter  contribution 


<o 


pQ> 


4  irk” 


■»pQ 


n*a 


P„(n^,n^|n)p(nJ  |xS(v*n) 


s  /"  _  —  \  1 2 


n-*n 


|  x8(v*n  ) | 2  <aPQ  > 


(2.24) 


in  which  x  (vn  )  =  1  is  the  characteristic  function  for  the  small  scale 
s 

PQ 

surface  height  roughness.h  ,  and  <a  >  is  the  physical  optics  scattering 

8  00 

cross  section  (Bahar  1981a, b) .  The  second  term  which  accounts  for  Bragg 
scattering  is  given  by 


<o 


W>. 


irk 


|DPV|2 

y — W(v-,v-)P„(n  ,n  |n)p(h  ,h  )dh  dh 


n*a 


X  Z  X  z 


(2.25) 


in  which  W(v-,v-)  is  the  small  scale  surface  height  spectral  density 
X  z 

(Barrick  1970,  Ishimaru  1978). 


W(v-,v-)  •  4-  I  <h  (x,z)h' (x' ,z ')>  exp(iv-x,  +  iv_z,)dx 
x*  x  2  a  *  s  *  r  x  d  x  d  < 

TT  y 


d  dzd  (2.26) 


In  (2.26)  it  is  assumed  that  the  small  scale  surface  autocorrelation 

function  depends  on  distances  [x^n^  +  zdn2^  measured  along  the  large 

scale  surface  y-h„(x,z)  ■  0,  and  h  ,  the  small  scale  surface  height  is 
.  **  s 

measured  penpendicular  to  the  large  scale  surface.  Furthermore 

;  -  +  v;»2  +  (2.27) 

is  the  expansion  of  the  vector  v  in  the  local  coordinate  system  (at 
any  point  on  the  large  scale  surface)  associated  with  the  unit  vectors 
np  ^2,  and  n^  (Bahar  1982a, b)  (see  Fig.  2.2). 


(n  x  a  )/|n  x  a 
z  z 


n2  ■  n,  n^  ■  tij  x  n 


(2.28) 


The  scattering  cross  section  (2.25)  is  in  complete  agreement  with  pertur¬ 
bation  theory  (Rice  1951,  Barrick  1970).  Thus<oP^>2  can  be  regarded  as 
an  average  (over  the  distribution  of  slopes  of  the  large  scale  surface) 
of  the  scattered  power  from  patches  of  slightly  rough  surfaces  that  ride 
the  large  scale  surface  (Wright  1966,  1968,  Valenzuela  1968).  On 
expressing  the  unit  vector  n  in  terms  of  the  slope  angles  and  6  in 
and  perpendicular  to  the  plane  of  incidence,  the  scattering  cross  section 
(2.25),  can  be  compared  with  earlier  solutions  that  are  "mostly  based 
on  physical  considerations"  (Valenzuela  1968,  Valenzuela,  Liang  and 
Daley  1971,  Bahar  1981c).  The  expression  (2.25)  is  also  in  agreement 
with  Brown's  solution  (that  is  expressed  in  terms  of  a  two-dimensional 

convolution  of  transforms)  provided  that  the  mean  squares  of  the  large 
2  2 

scale  slopes  o__  ■  o  are  very  small.  The  difference  between  Brown's 
X  z 


solution  and  the  full  wave  solution  arises  because  Brown  (1978)  (on 


incidence  and  scatter  and  local  coordinate  system  with  unit  vectors 


using  a  combination  of  Burrows'  perturbation  theory  (1967)  and  physical 
optics  (Beckmann  1968)  assumes  that  the  surface  height  autocorrelation 
function  for  the  small  scale  surface  is  dependent  on  distances  in 
the  reference  (mean)  plane  rather  than  distances  along  the  large  scale 
surface,  as  assumed  in  this  work  and  implicitly  assumed  by  Valenzuela 
(1968).  Furthermore,  it  should  be  pointed  out  that  in  Burrows'  pertur¬ 
bation  theory  the  small  scale  surface  height  is  the  distance  from  the 
unperturbed  (filtered)  surface  to  the  perturbed  surface  measured  along  a 
line  perpendicular  to  the  unperturbed  surface  rather  than  perpendicular 
to  the  reference  surface  as  assumed  by  Brown. 

Since  the  full  wave  solution  (2.23)  accounts  for  both  Bragg  scat¬ 
tering  as  well  as  specular  point  scattering  in  a  self-consistent  manner 
without  introducing  a  combination  of  perturbation  and  physical  optics 
theories,  it  is  not  necessary  to  decompose  the  surface  h^  into  large 

and  small  scale  surface  heights  h9  and  h  in  order  to  evaluate  the  total 

S 

PQ 

scattering  cross  section  <a  >.  Moreover,  it  is  not  necessary  to 
restrict  the  application  of  the  full  wave  approach  only  to  surfaces  that 
can  be  regarded  as  a  small  scale  surface  perturbation  superimposed  on  a 
large  scale  surface  roughness  with  very  large  radii  of  curvature  compared  to 
the  wavelength.  Thus  in  Section  C.4)the  full  wave  approach  is  applied  to 
more  general  composite  rough  surface  models  that  may  frequently  fail  to 
satisfy  these  restrictions. 

2.4  Application  to  Surfaces  That  Do  Not  Satisfy  Perturbation  and 


composite  surface  h(x,z),  the  total  scattering  cross  section  can  be 
evaluated  using  the  full  wave  solution  (2.18).  Thus  using  the  full 


wave  approach  it  is  not  necessary  to  specify  the  wave  number  kj  at 

which  spectral  splitting  between  the  two  scale  surface  heights  is 

assumed  to  occur.  On  the  other  hand,  when  different  theories  such 

as  physical  optics  and  perturbation  theories  are  combined  to  analyze 

composite  rough  surfaces,  it  is  necessary  to  spectrally  decompose 

(filter)  the  rough  surface.  For  instance,  in  Brown's  work  (1978)  the 

specification  of  k^  is  assumed  to  be  based  upon  the  characteristics 

2  2 

of  the  small  scale  structure  (k  <h  >  <<  1).  However,  in  the  work  by 

os  J 

Tyler  (1976),  the  specification  of  kj  is  assumed  to  be  based  on  the 
characteristics  of  the  large  scale  surface.  Thus,  in  order  to  apply 
physical  optics  theory  to  the  large  scale  surface  roughness,  Tyler 
imposes  the  condition 


2  *9  2 

in  which  P^r^l  <  is  t*ie  probability  that  Ir^l  (the  absolute 

value  of  the  product  of  the  principal  radii  of  curvature  of  the  large 

"2  '2 

scale  surface)  is  less  than  r^»  and  the  critical  value  for  r^  is 

assumed  to  be  k  ^  (Tyler  1976).  If  both  conditions  (2.22b)  (k^<h^>  <<  1) 
o  os 

2  2 

as  well  as  (2.29)  (P( | r 12 |  <  r^2)  <<  1)  are  satisfied  simultaneously, 

the  scattering  cross  section  is  expressed  as  the  sum  of  the  specular 
point  scattering  cross  section  and  the  Bragg  scattering  cross  section 
(2.23),  since  surfaces  h  and  h  individually  satisfy  the  limitations 
imposed  by  perturbation  and  physical  optical  theories  respectively. 
However,  if  both  conditions  are  not  satisfied  simultaneously,  the 


physical  optics  model  or  the  perturbation  model  or  even  a  perturbed- 
physical  optics  model  cannot  be  assumed  in  general. 

To  gain  more  physical  insight  and  to  facilitate  the  evaluation  of 
the  cross  section  (2.1),  it  is  assumed  here  that  the  rough  surface 
height  h  is  decomposed  into  two  surfaces  h^,  and  h^  (see  Fig.  2.3).  Thus 
r  -  rp  (x,tyz)  +  lyi  (2.30) 

in  which,  following  the  filtering  scheme  proposed  by  Tyler  (1976),  the 

2  *2  2  2 

surface  h^,  satisfies  the  conditions  P  ( I r  ^2 1  <  rl2^  <<:  *  an<*  >:>  * 

However,  it  is  assumed  that  the  remainder  term  h^  does  not  satisfy  the 

2-2 

perturbation  condition  (kQ<h^>  <<  1).  Assuming  that  hp  and  h^  are 
statistically  independent  random  functions 


<exp£ivy(h-h,)]>  -►  exp£ivy(h^xd+h^zd)]X2(v-,-v-) 


(2.31a) 


Thus  the  full  wave  solution  (2.18)  can  be  written  as 
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SPQexp[ivy(h^xd+h^zd)](|xR|2  +XR"|xR|2)  (2.31b) 


exP['ivxxd  +  ivzzd]dxd  dzd> 


Ri2 


in  which  for  convenience  the  term  |x  |  is  added  and  subtracted  and  x 
R 


and  X2  are  characteristic  and  joint  characteristic  function  for  the 
surface  height  hj^  respectively. 

(2.31c) 


X2(vy,_Vy>  "  <exPtiv-(hR-h^)]> 


F  0  2  2 

In  (2.31b)  it  is  assumed  that  |x  |  <<  1  since  kQ<hp>  >>  1.  Thus,  fol- 

lowing  the  analytical  procedures  used  in  deriving  (2.23)  (Bahar  1981c), 


it  can  be  shown  that 


rd  =  DISTANCE  MEASURED  ALONG  LARGE  SCALE  SURFACE 


47 


in  which 


k2  f 


SPQexpfiVy(h^xd+h^d)]|xR|2exPrivxxd+iv2zd]dxddzd> 


lxR^-no)  )2  <oPQ> 


(2.32b) 


and  <o  >  is  the  specular  point  scattering  cross  section  for  the 


filtered  surface  h^.  Since  in  this  case  it  is  assumed  that  k2<h^> 
not  much  smaller  than  unity,  the  factor  in  (2.32),  |xR|2,  can  be  signi¬ 


ficantly  different  from  unity.  On  deriving  the  expressions  for  <or^> 


PO  2  2 

(2.24),  and  <a  >  (2.32b),  it  is  assumed  that  v  <h  >  >>  1  and  deep 

r  y  r 


phase  modulation  occurs  for  all  the  roughness  scales  included  in  h  . 

F 


Furthermore,  the  scattering  cross  section  associated  with  the  sur¬ 


face  h^  (that  rides  on  the  filtered  surface)  is 
nn  ko  f  I  °PQ  1 2po  (nf  .n^n)  _  _ 


CxJCVy.-  v">  “  ix  (v-)|  ]  (2.32c) 


• expE iv_  xd  +  iv-zd]dxd  dzd> 


in  which  it  is  assumed  that  (2.19)  is  valid,  v_,  v-  and  v-  are  the  com- 

x  y  z 


ponents  of  v  in  the  local  coordinate  system  associated  with  the  filtered 


-2  -2  k 

surface  h^  (2,27),  and  (xd  +  zd'  *s  distance  measured  along  the 


filtered  surface  y-h_(x,z)  -  0  (see  Fig.  2.3). It  can  be  readily  shown 


that  if  kQ<h^>  <<  1  (2.32c)  reduces  to  (2.25)  with  h^  *  hg,  since  in  this 


x2<vy--y  -ixR<y  i2  ■  <E*  y 

If  the  surface  height  h^  is  normally  distributed  and 


(2.33) 
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ko<H^>  is  not  much  smaller  than  unity 


X^Vp.-v-)  -  lxR(Vp)  1 2_exp(~Vy<h^>) [exp(v^<hRh^>)~  l].  (2.34 


In  this  case  (2.34)  is  expanded  as  follows 

u  tj  o  7—7  00 

x2(vy»~vy}  -lx  (v-)|  =  exp(-v-<i^>)E^  - - 

and  <oPQ>D  (3.22c)  can  be  expressed  as 
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(2.35) 


<oPQ>. 


(2.36) 


where 
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(2.37) 
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(2.38) 


In  (2.38)  the  symbol® denotes  the  two  dimensional  convolution  of  W  , 
W, (v-,v-)  W(v_,v_) 

with  W.  Since  — — —  ■  - 7 -  is  the  two  dimensional  Fourier 

4  4 

transform  of  the  surface  height  autocorrelation  function  <hRh^>,  the 


term  <crw>R^  accounts  for  first-order  Bragg  scattering.  However,  since 

2  -o  2  —2 

4k  <1£>  is  not  much  smaller  than  unity  the  factor  exp(-v_  h^  )  appear- 

O  K  y  K 

PO 

ing  in  <o  X>R^  is  significantly  smaller  than  unity  for  backscatter  near 
normal  incidence  and  approaches  unity  for  backscatter  near  grazing 
incidence.  This  factor  does  not  appear  in  the  solution  derived  on  the 
basis  of  perturbation  theory  (2.25). 

When  the  autocorrelation  function  <hRH^>  can  be  expressed  analytic¬ 
ally  in  closed  form,  the  two  dimensional  Fourier  transforms  in  (2.38) 
can  be  integrated  directly.  For  example,  when  the  Gaussian  form  is 


used. 


<W  =e*p<(-*d  ‘  *d)/l2) 


(2.39) 


closed  form  expression  can  be  derived  for  W  (v  ,v_)  which  are  also 

ID  X  Z 

Gaussian  (Beckmann  and  Spizzichino,  1963).  However,  since  in  practical 

problems  the  surface  height  is  usually  characterized  explicitly  by 

its  spectral  density  function  W(v-,v-),  and  since  this  function  is  rarely 

x  z 

Gaussian  for  natural  surfaces,  the  Fourier  transforms  in  (2.38)  are 
evaluated  through  the  repeated  use  of  the  two  dimensional  convolution 
theorem. 

If  the  mean  square  slope  of  the  surface  height  is  very  small  com¬ 
pared  to  unity  (n  -►  ay) 


(2.40) 
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The  expression  <a  H>j£  can  readily  evaluated  for  backscatter 

near  normal  incidence.  Thus  assuming  a  perfectly  conducting  surface, 

,  -f  -i 
for  n  ■  -n  -  a 


<aP\2  -*•  exp(-4k^<h^>)W2(0,0)  (2.41) 

Assuming,  for  example,  as  isotropic  surface  height  spectral  density 


of  the  form 


kj  <  k  <  k2 


k  <  k^  and  k  >  k2 


(2.42) 


2*  k2  2 
W2(0,0)  m  f  T  W2 (k)kd  kd<|>  =  J-J 

°  kl  U  k2  . 

The  mean  square  of  the  surface  height  h^  is 
2ir  k„ 
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Thus  for  k2  >>  k2 

{•A2  4<iV  3 

V°-°>  -  R - ST- 


and  (2.41)  reduces  to 


<oPP>R2  "  k  (^o^>)3exp(-4k2<^» 


(2.43) 


(2.44) 


(2.45) 


(2.46) 


As  in  (2.25),  the  contribution  (2.32c)  to  the  scattering  cross 

section  is  due  to  scattering  by  the  surface  that  rides  on  the  filtered 

2  -2 

surface  hp.  In  general  however,  if  kQ<h^>  Is  not  much  smaller  than  unity 


PO 

the  contribution  to  the  scattering  cross  section  <a  x>  ,  (2.32c) 

K 

cannot  be  derived  on  the  basis  of  perturbation  theory  and  the  factor 

|X  I  in  (2.32b)  cannot  be  replaced  by  unity.  Thus  a  combination  of 

perturbation  theory  and  physical  optics  theory  cannot  be  applied  to 

R  2 

this  problem.  Physically,  the  coefficient  |x  |  in  (2.32b)  accounts 

2  -2 

for  the  fact  that  when  kQ<h^>  is  not  much  smaller  than  unity,  the 

that  ride  on  the  filtered  surface  hp  could 
significantly  reduce  the  contribution  of  the  "specular"  scattering 

PO 

cross  section,  <a  >^,  associated  with  the  filtered  surface  h^.  For 

backscatter,  the  perturbation  term  associated  with  Bragg  scatter 
PQ 

<o  (2.25)  is  much  smaller  than  the  term  associated  with  specular 

scatter  <qP^>  (2.24)  near  normal  incidence  and  <0PP>, »  <aPP>  for 

O'  1  o 

near  grazing  Incidence  provided  that  (2.22)  and  the  radius  of  curva¬ 
ture  criteria  are  satisfied.  However,  for  the  more  general  case  treated 

PP 

in  this  section  a  similar  relationship  between  <o  >  (2.32b)  and 

F 

PP 

«j  may  not  exist.  Thus,  while  it  is  certainly  not  necessary  to 
decompose  (filter)  the  rough  surface  height  h  into  surfaces  hp  and 
(2.30)  in  order  to  derive  the  scattering  cross  sections  when  the  full 
wave  approach  is  used,  such  a  decomposition  does  provide  additional 
physical  insight.  The  results  of  the  analysis  carried  out  in  this 
section  illustrate  how  the  specular  scattering  and  Bragg  scattering  com¬ 
ponents  of  the  total  scattering  cross  section  begin  to  blend  with  each 

2  -2 

other  as  the  value  of  kQ<h^>  increases  and  the  composite  surface  can 
no  longer  be  regarded  as  a  perturbed-physical  optics  (Kirchhoff)  surface. 
Furthermore,  the  decomposition  (filtering)  of  the  surface  height  assists 
in  making  the  result  (2.1)  easier  to  compute. 


surface  irregularities 


It  can  also  be  shown  that  If  Instead  of  (2.30)  the  criteria 
for  decomposing  the  composite  surface  is  based  only  on  the  character¬ 
istics  of  the  small  scale  surface  height  h  ,  (2.22b)  such  that 

s 

h  -  hs  +  \  ».«> 
and  if  the  radii  of  curvature  associated  with  the  remaining  surface 
are  not  large  compared  to  wavelength  (such  that  the  physical  optics 
theory  cannot  be  applied  to  it)  the  full  wave  approach  can  still  be  used. 
2.5  Concluding  Remarks 

It  is  shown  that  since  the  full  wave  theory  accounts  for  both  Bragg 
scattering  as  well  as  specular  point  scattering  in  a  self-consistent 
manner,  in  order  to  evaluate  the  scattering  cross  sections  it  is  not 
necessary  to  use  a  combination  of  perturbation  and  physical  optics  theo¬ 
ries.  Furthermore,  it  is  shown  that  the  general  formulas  derived  in 

Section  (2.3)  can  be  applied  to  surfaces  that  do  not  necessarily  satisfy 
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the  perturbation  theory  restrictions  kQ<h^>  <<  1  and/or  the  physical 

optics  theory  restrictions  on  the  radii  of  curvature  (2.29).  For  the 
general  case  considered  in  Section  (2.4),  in  which  a  perturbed-physical 
optics  approach  cannot  be  used,  it  is  illustrated  how  Bragg  scattering 
and  specular  point  scattering  begin  to  blend  with  each  other.  In  this 
case,  decomposition  (filtering)  not  only  enhances  one's  physical  insight, 
but  also  facilitates  the  numerical  evaluation  of  the  scattering  cross 
sections  (2.1).  These  numerical  evaluations  need  to  be  pursued  with 
special  emphasis  on  the  specification  of  (where  spectral  splitting 
is  assumed  to  occur) . 


3.0  COMPUTATIONS  OF  SCATTERING  CROSS  SECTIONS  FOR  COMPOSITE  SURFACES 
AND  THE  SPECIFICATION  OF  THE  WAVENUMBER  WHERE  SPECTRAL  SPLITTING 
OCCURS 


3.1  Background 


The  scattering  cross  sections  for  composite  random  rough  surfaces 
are  evaluated  using  the  full  wave  approach.  They  are  compared  with 
earlier  solutions  based  on  a  combination  of  perturbation  theory  which 
accounts  for  Bragg  scattering  and  physical  optics  which  accounts  for 
specular  point  theory.  The  full  wave  solutions  which  account  for  both 
Bragg  scattering  and  specular  point  scattering  in  a  self-consistent 
manner  are  expressed  as  a  weighted  sum  of  two  cross  sections.  The  first 
is  associated  with  a  filtered  surface,  consisting  of  the  larger  scale 
spectral  components,  and  the  second  is  associated  with  the  surface 
consisting  of  the  smaller  scale  spectral  components.  The  specification 
of  the  surface  wavenumber  that  separates  the  surface  with  the  larger 
spectral  components  from  the  surface  with  the  smaller  spectral  compon¬ 
ents  is  dealt  with  in  detail.  Since  the  full  wave  approach  is  not 
restricted  by  the  limitations  of  perturbation  theory,  it  is  possible 
to  examine  the  sensitivity  of  the  computed  values  for  the  backscatter 
cross  sections  to  large  variations  in  the  value  of  the  wavenumber  where 
spectral  splitting  is  assumed  to  occur. 

3.2  Discussion 


In  this  section,  the  backscatter  cross  sections  for  composite 


models  of  rough  surfaces  are  evaluated  using  the  full  wave  solutions 
to  the  problem  (Bahar  and  Barrick  1982).  The  like-polarized  backscatter 
cross  sections  for  both  vertically  and  horizontally  polarized  waves  as 


well  as  the  cross-polarized  backscatter  cross  sections  are  evaluated 
as  functions  of  the  angle  of  incidence.  In  order  to  compare  the  full 
wave  results  with  earlier  solutions  for  the  backscatter  cross  sections, 
the  full  wave  solutions  are  expressed  in  terms  of  a  weighted  sum  of 
scattering  cross  sections.  The  first  is  associated  with  a  filtered 
surface  h^,  consisting  of  the  larger  scale  spectral  components  (k  <  kj)  . 
The  second  is  associated  with  a  surface  h^  consisting  of  the  smaller 
scale  spectral  components  (k  >  kj)  (see  Section  2). 

In  an  attempt  to  draw  more  definite  conclusions  regarding  the 


choice  of  k^  (the  wavenumber  where  spectral  splitting  is  assumed  to 

occur)  between  the  surfaces  h^,  and  h^,  the  wavenumber  kj  is  varied  over 

a  wide  range  of  values.  The  wavenumber  k^  is  related  to  the  parameter 
2  2 

8  ■  4kQ  <h^> ,  where  kQ  is  the  wavenumber  for  the  electromagnetic  wave 
2 

and  <h^>  is  the  mean  square  of  the  surface  height  h^.  Thus,  on  applying 

a  perturbed-physical  optics  approach  to  rough  surface  scattering  (Brown 

1978,  1980),  the  wavenumber  k,  was  chosen  on  the  basis  of  the  charac- 

d 

teristics  of  the  surface  consisting  of  the  small  scale  spectral  compon¬ 


ents.  However,  in  the  approaches  of  Hagfors  (1966)  and  Tyler  (1976), 


the  specification  of  the  wavenumber  kj  is  assumed  to  be  based  on  the 
characteristics  of  the  filtered  surface  h  .  In  their  approaches,  however, 
they  ignore  the  effects  of  the  surface  (consisting  of  the  smaller 
spectral  components,  k  >  kj)  since  their  results  were  meant  to  explain 
backscatter  from  lunar  and  planetary  surfaces  at  near  normal  incidence. 


It  should  be  pointed  out  that  even  if  the  spectral  components  of 
the  filtered  surface  density  satisfy  the  radii  of  curvature  criteria 
(associated  with  the  Kirchhoff  approximations  of  the  surface  fields). 
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they  may  not  necessarily  satisfy  the  condition  for  deep  phase  modulation 

Implicit  in  the  evaluation  of  the  specular  point  result  (Barrick  1970). 

In  Section  (3.3)  of  this  report  the  principal  full  wave  expressions 

for  the  normalized  scattering  cross  sections  are  summarized.  Since  the 

parameter  g  need  not  be  restricted  by  considerations  of  perturbation 

theory,  it  is  assumed  that  the  filtered  surface  h,  satisfied  the  radii 

F 

of  curvature  criteria  (associated  with  the  Rirchhoff  approximations  for 
the  surface  fields)  as  well  as  the  conditions  for  deep  phase  modulation. 

In  Section  (3.4)  an  extensive  set  of  numerical  data  is  presented 
for  the  backscatter  cross  sections  associated  with  the  like-polarized 
and  the  cross-polarized  waves.  It  is  shown  that  while,  as  expected,  the 
cross  sections  associated  with  the  individual  surfaces  h^  and  h^  criti¬ 
cally  depend  on  the  choice  of  3  (and  therefore  kj) ,  the  total  back¬ 
scatter  cross  section  remains  practically  insensitive  to  6  for  3  _>  1.0. 
However,  there  are  small,  though  perhaps  significant,  differences  in  the 
values  for  the  total  cross  sections  where  3  is  increased  from  0.1  to  1.0. 
It  should  be  pointed  out  that  these  differences  (as  3  varies  from  0.1 
to  3  ■  1.0)  are  significantly  smaller  than  those  predicted  on  the  basis 
of  Brown's  analysis.  The  reason  why  the  full  wave  results  for  the  total 
backscatter  cross  sections  merge  for  values  of  3  M..0  is  related  to  the 
condition  for  deep  phase  modulation.  Some  of  the  key  observations  of 
Brown,  in  his  contribution  to  rough  surface  scattering,  (1978),  are 
summarized  in  Section  (3.4)  to  emphasize  the  problems  related  to  the 
specification  of  kj.  However,  the  interested  researcher  in  this  field 
should  familiarize  himself  with  his  work  as  well  as  the  pioneering 


■■.y. 


Va 


Ji 


•/VV  v  v" 


contributions  of  Hagfors  (1966)  and  Tyler  (1976). 


The  details  given  in  the  illustrative  examples  of  Section  (3.4)  are 
presented  primarily  for  tutorial  purposes  and  to  vividly  establish  the 
criteria  for  specifying  kj.  Thus,  the  engineer  need  not  consider  the 
evaluation  of  the  cross  section  for  6  >  1.0  or  concern  himself  with  the 
individual  terms  that  add  up  to  determine  the  total  cross  sections. 
Nevertheless  it  does  enhance  the  engineer's  physical  insight  in  dealing 
with  problems  of  rough  surface  scattering.  For  example,  the  full  wave 
approach  can  explain  why  the  frequency  dependence  of  the  backscatter 
cross  sections  changes  as  one  changes  the  angle  of  incidence  or  when  one 
varies  the  frequency  of  the  electromagnetic  wave. 

3.3  Formulation  of  the  Problem 

In  this  section  the  principal  expressions  for  the  normalized  scat¬ 
tering  cross  sections  per  unit  area  are  summarized  and  the  full  wave  solu 
tions  are  compared  with  earlier  solutions  based  on  perturbation  and 
physical  optics  theories. 

The  expressions  for  the  normalized  scattering  cross  sections  per 
unit  area  (Ishimaru  1978),  based  on  the  full  wave  solutions  for  the 
incoherent  scattered  radiation  fields,  are  given  by  (Bahar  1981a) 


<«PQ>  -  £  /  [<sPQ 


exp[ivy(h-h')]> 


DP<^P2(nf  .n^n)  .1 

'  I  - 2— -  >  )  1  J 

n*a  3 


•exp[iv  x,  +  iv  z.ldx.dz, 
r  x  d  z  d  d  d 


(3.1) 


(3.2) 


in  which 

r ,  ■  (x-x')a  +  (z-z')a  ”  x.a  +  z.a 

d  x  '  z  dx  dz 

is  the  radius  vector  between  two  points  in  the  reference  plane  (x,z) 
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(see  Fig.  3.1).  The  vector  v  in  the  reference  coordinate  system  (x,y,z) 


v  -  k  (n^-n1)  “  v  a  +va  +  va 
o  xx  y  y  2  z 


(3.3) 


where  kQ  is  the  free  space  wavenumber  for  the  electromagnetic  wave  and 
n^  and  are  unit  vectors  in  the  direction  of  the  incident  and  scattered 
wave  normals  respectively.  An  exp  (iut)  time  dependence  is  assumed 
throughout  this  work.  The  symbol  <  >  denotes  the  statistical  average 
and  it  can  be  shown  that  (Bahar  1981b) 

f  PQ  -  ^ 

<SPQ  exp[iv  (h-h')]>-  <  1  ~  ^  exp[iv  (h-h')] 
y  )  n*a  y 

J  V 


p(hx,hz)P2(nf,5i|n)dhxdh2  > 


(3.4) 


where  n  is  the  unit  vector  normal  to  the  rough  surface  (see  Fig.  3.1). 


f(x,y,z)  -  y  -  h(x,z)  -  0 


Vf  -  n  Vf  -  V 


|y-h(x ,  z)| 


(-h  a  +  a  -ha) 
xx  y  z  z 


in  which  the  components  of  the  gradient  of  h(x,z) 

h  ■  3h/3x  ,  h  ■  3h/3z 

x  z 


(3.5a) 


(3.5b) 


(3.5c) 


are  random  variables  and  p(h  ,h  )  is  the  distribution  function  for  the 

X  2 

slopes  h  and  h  .  The  expression  for  <oP<^>  (3.1)  accounts  for  shadowing 


P  (n  ,n  |n)  -  P,(n  ,51|H  )S(n  •5)S(-n1*H) 


(3.6) 


in  which  P2(n^,n^|n)  is  the  probability  that  a  point  in  the  rough  surface 
is  both  Illuminated  and  visible  given  the  value  of  the  slopes  at  the  point 
(Smith  1967;  Sancer  1969)  and  P2(n^,n*|n  )  is  its  value  at  the  specular 


points  where  the  unit  vector  n  is  given  by 


5  *  '  T?r?iT  ■  '?/v  (3-7) 

jn  -n  | 

The  arguments  of  the  unit  step  functions  S(-ni*n)  and  S(nf*n)  vanish  at 

points  of  the  rough  surface  where  the  incident  and  scattered  waves  are 

tangent  to  the  surface.  Thus  S(-n^*n  )  ■  1  and  S(ii^*n  )  *  1 .  The  charac- 

s  s 

teristic  function  and  the  joint  characteristic  function  for  the  surface 
height  h  are  respectively, 

X(vy)  -  <exp(ivyh)>  (3.8) 

and 

x2 (vy > ”vy)  ■  <exp  iVy(h-h’)>  (3.9) 

For  rough  surface  heights  with  Gaussian  distributions  (assumed  in  this  work) 

X(vy)  -  exp(-Vy<h2>/2)  (3.10) 
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in  which  <h  >  -  a  is  the  mean  square  surface  height.  Furthermore, 

X2^vy)  “  exp[-v2(<h2>  -  <h  h * >) ]  (3.11) 

where <h(x,z)h' (x’ ,z’)>  is  the  surface  height  autocorrelation  function.  The 
coefficients  D  x  (Bahar  1981a)  depend  explicitly  upon  the  polarization  of 
the  incident  wave  (second  superscript  Q*V  -  vertical,  Q-H  -  horizontal) 
and  the  polarization  of  the  scattered  wave  (first  superscript  (P*»V,H), 
the  direction  of  the  incident  and  scattered  wave  normals  n*  and  n^  re¬ 
spectively,  the  complex  permittivity  and  permeability  of  the  medium  of 
propagation  c  and  y  respectively  and  the  unit  vector  n  normal  to  the  rough 
surface.  In  the  above  expression  for  <o  >,  (3.1),  it  has  been  assumed 
that  the  surface  height  h  and  slopes  (^x»^z)  are  statistically  independent 


a  condition  that  holds  for  Gaussian  surfaces  at  each  point)  and  that  for 
distances  |r^|  less  than  the  correlation  distances  £  (where  the  surface 
height  autocorrelation  function  C(£)  -  C(o)/e) 

n(h  ,h  )  =  n' (h' ,h')  (3.12) 

X  Z  X  z 

If  the  statistics  for  the  entire  rough  surface  h(x,z)  are  known,  the 
total  scattering  cross  sections  can  be  evaluated  using  the  full  wave  solu¬ 
tion  (3.1).  However,  when  a  combination  of  perturbation  theory  (Rice 
1951;  Barrick  1970)  and  physical  optics  theory  (Beckmann  1968)  is  applied 
to  the  problem  of  rough  surface  scattering,  a  two  scale  model  of  the 
surface  is  used  and  the  surface  is  decomposed  into  a  filtered  surface 
height  (consisting  of  the  large  scale  spectral  components  of  the  surface 

height)  and  a  small  scale  surface  height  h  that  is  superimposed  (rides  on) 

s 

the  large  scale  filtered  surfaces  (Wright  1968,  Valenzuela  1968;  hg 
consists  of  the  smaller  scale  spectral  components  of  the  surface  height) . 

To  this  end  it  is  necessary  to  specify  the  wavenumber  kj  at  which  spectral 
splitting  is  assumed  to  occur.  For  instance.  Brown  (1978)  who  uses  a  com¬ 
bination  of  Burrows'  perturbation  theory  (1967)  and  physical  optics 

(Beckmann  1968),  bases  the  specification  of  kj  upon  the  characteristics  of 

2  2 

the  small  scale  structure  (k  <h  >  <<  1).  However,  in  the  works  by 

os 

Hagfors  (1966)  and  Tyler  (1976),  the  specification  of  kj  is  assumed  to  be 

based  upon  the  characteristics  of  the  large  scale  (filtered)  surface  h  . 

F 

Thus,  in  order  to  justify  the  application  of  physical  optics  theory  to  the 
large  scale  surface  roughness,  Tyler  (1976)  imposes  the  condition 


in  which 


pc  Ir12 1  <  ^12)  <<:  1 
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P(  | r ^ 2 1  <ri2^  is  t*ie  Probability  that 


(3.13) 


(the  absolute  value 
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of  the  product  of  the  principal  radii  of  curvature  of  the  large  scale 

„2  „2 

surface)  is  less  than  r^*  an<*  critical  value  for  r^  Is  assumed  to  be 
kQ'2 (Tyler  1976). 

Since  the  full  wave  approach  accounts  for  both  specular  point  scatter¬ 
ing  as  well  as  Bragg  scattering  in  a  self-consistent  manner,  it  is  not 
necessary  to  filter  (decompose)  the  rough  surface  to  evaluate  the  scatter¬ 
ing  cross  sections.  However,  filtering  the  composite  surface  enhances 
one's  physical  insight  as  to  the  validity  (or  lack  thereof),  of  the 
perturbed -physical  optics  approach  to  the  scattering  problem  and  also 
facilitates  the  numerical  evaluation  of  the  cross  section.  Thus  it  is 
assumed  here  that  the  rough  surface  height  h(x,z)  is  decomposed  into  two 
surfaces  such  that  the  position  vector  to  a  point  on  the  rough  surface  is 
rs  “  r^x.hp.z)  +  n  (3.14) 

In  order  to  apply  the  physical  optics-specular  scattering  approximation 
to  the  filtered  surface  h^,  (consisting  of  the  larger  scale  spectral  com¬ 
ponents)  assume  that  h^,  satisfies  the  appropriate  radii  of  curvature  cri¬ 
teria  (for  example  (3.13)).  In  addition,  assume  that  deep  phase  modulation 
occurs.  Thus,  the  distances  from  the  transmitter  and  receiver  to  the  indi¬ 
vidual  specular  points  (which  are  random  variables)  are  such  that  the 
contributions  from  the  individual  specular  points  are  distributed  uniform¬ 
ly  in  phase  from  -ir  to  it..  One  should  also  note  that  no  matter  how  the  com¬ 
posite  surface  is  filtered,  the  physical  optics  approach  (based  on  the 
Kirchhoff  approximations  for  the  surface  fields)  is  not  valid  if  for  a 
given  n*  and  n^  specular  points  do  not  exist  on  the  rough  surface 


r.  s  v*. 
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is  specified  by  the  desired  characteristics  of  the  filtered  surface  h^, it 

2  -2 

is  not  assumed  here  that  kQ<hR>  is  much  smaller  than  unity.  Thus  per¬ 
turbation  theory  (Rice  1951;  Barrick  1970),  cannot  be  applied  to  the 
remaining  part  of  the  surface  height  h^.  For  truly  random,  rough  natural 

surfaces  such  as  the  sea  surface,  it  is  assumed  that  h  and  h  are  statis- 

t  K 

tically  random  functions.  If  in  addition,  it  is  assumed  that  Vh  -  Vh  , 

r 

the  full  wave  solution  (3.1)  can  be  expressed  as  a  weighted  sum  of  the 
individual  cross  sections  for  the  surfaces  h^  and  hR  respectively  (Bahar 


1981b). 


PQ  PQ  ,  PQ 

<a  *»  <a  x>  +  <a 

F  R 


(3.15) 


The  first  term  <o  can  be  shown  to  be  given  by 
<aPQ>  -  |xR(v*n)|  <oPQ>  , 


(3.16) 


R  — •  —  R 

in  which  x  (vn  )  *  x  (v)  is  the  characteristic  function  for  the  surface 
8 

PQ 

hR  and  <o^  >  is  the  specular  point  scattering  cross  section  for  the  fil¬ 
tered  surface  h  .  The  factor  XR(v)  that  multiplies  <oP<^>  accounts  for  the 
degradation  of  the  specular  points  contributions  due  to  the  superimposed 
surface  h^  (Bahar  1981a).  It  can  be  shown  that 


>PQ>  -  — 2. 
00  v2 

y 


P2(nf  ,ni|n)P(n) 

'n-a  1  n-+ii 

y  J  s 


(3.17) 


The  second  term  <a  is  the  scattering  cross  section  for  the  surface 
that  rides  on  the  filtered  surface  hy.  It  can  be  expressed  as  follows 
(Bahar  and  Barrick  1982) 


PQ  " 

<0  \  "  I 


(3.18) 


<0 


pQ> 


4irk 


m 


f  l»  I  |!>  ,  2  r2  i  pj-  Wm(v-,v-) 


n*a 


exp(-v-  <hR>) 


p(h  ,h  )dh  dh 
r  X9  Z  X  z 


(3.19) 


in  which  v-,  v-  and  v-  are  the  components  of  v  (3.3)  in  the  local 


coordinate  system  (at  each  point  on  the  large  scale  surface)  associated 
with  the  unit  vectors  n^ ,  n 2,  and  n^  (see  Fig.  3.2).  Thus  v 
in  (3.3)  is  also  expressed  as  (Bahar  1982) 


v  “  v-n.  +  v-n0  +  v-n, 
x  1  y  2  z  3 


(3.20) 


where 


■  (n  x  az)/|n  x  aj  ,  n2  -  n  ,  n3  -  irij  x  n 


(3.21) 


.2m 


The  function  W  (v-,v-)/2  is  the  two  dimensional  Fourier  transform  of 
m  x  z 


W  (v-,v-)  . 

m  x  z  1 


,2m 


(2ir) 

_ l 

_2m 


/  «W>* 


exp(iv-xd  +  iv-zd)dxddid 


- — f  W  . (v-' ,v-' 
m  I  m-1  x  z 


)W.  (v-  -  v- ' ,v-  -  v-')  dv-'dv-' 
l  x  x  z  z  x  z 


-  W  ,(v-,V-)(jT)  W.(v-,v-) 

22m  m-1  x  z  v—/  l  x  z 


(3.22) 


In  (3.22)  |x,a  +  z  a  ]  is  the  distance  measured  along  the  large  scale 
Q  x  a  z 

surface  and  the  symbol  (xj  denotes  the  two  dimensional  convolution  of  W 

' — '  m- 

with  W, .  Since  W.(v-,v-)/4  -  W(v-,v-)/4,  is  the  two  dimensional  Fourier 
1  1  x  z  x  z 

transform  of  the  surface  height  autocorrelation  function  <hRh^>,  it  is 
equal  to  the  spectral  density  for  the  surface  height  hR.  Thus,  the  first 
term  in  (3.18),  <o^>_.f 


accounts  for  first  order  Bragg  scattering. 


Incidence  and  scatter  and  local  coordinate  system  with  unit  vectors 
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However,  when  4k  <h£>  is  not  much  smaller  than  unity  the  factor 

2  -2  PQ 

exp(-v-<hR>)  appearing  in  the  expression  for  <o  X>R^  Is  significantly 

smaller  than  unity  for  backscatter  near  normal  incidence  and  approaches 


unity  for  backscatter  near  grazing  incidence.  This  factor  as  well  as  the 
R  PO 

coefficient  x  (v)  in  the  expression  for  <a  x>  (3.16)  do  not  appear  in 

r 

the  expressions  for  the  scattering  cross  sections  based  on  a  perturbed- 


physical  optics  solution.  Moreover,  the  arguments  of  W(v-,v-)  are  the 

x  z 

components  of  v  in  the  local  tangent  plane  and  not  in  the  reference  plane 
2  -2  PO  PO 

(x,z).  Thus  for  k  <h  >  <<  1,  <o  -  <o  x>  .  can  be  regarded  as  an  aver- 
O  R  R  R1 


age  (over  the  distribution  of  slopes  of  the  large  scale  surface)  of  the 


scattered  power  from  patches  of  slightly  rough  surfaces  that  ride  the 

large  scale  surface  (Wright  1966,  1968;  Valenzuela  1968).  The  scatter- 

PO 

ing  cross  section  <a  X>R1  can  be  compared  with  earlier  solutions  that 
are  "mostly  based  on  physical  considerations"  on  expressing  the  unit 
vector  n  in  terms  of  the  slope  angles  and  6  in  and  perpendicular  to  the 

plane  of  incidence  (Valenzuela  1968;  Valenzuela,  Laing  and  Daley  1971; 

2  -2  PQ 

Bahar  1982b).  For  kQ<h^>  <<  1  the  expression  for  <o  >R1  is  also  in 

agreement  with  Brown  (1978)  (as  corrected  in  Brown  1980) ,  provided  that  the 

mean  squares  of  the  large  scale  slopes  are  small.  This  difference  arises 

primarily  because  Brown  assumes  that  the  surface  height  autocorrelation 

function  for  the  small  scale  surface  is  dependent  on  distances  J r^ |  = 

2  2  k 

(x,  +  z,)  (3.2)  in  the  reference  (mean)  plane  rather  than  distances 

d  d 

—2  —2  h 

(x‘  +  z  S*  measured  along  the  large  scale  surface  as  assumed  in  this  work 
d  d 


and  implicitly  by  Valenzuela  (1968).  Moreover,  using  the  perturbation 


method  (Burrows  1967),  the  small  scale  surface  height  hR  is  measured 

normal  to  the  filtered  surface  h^  rather  than  normal  to  the  reference 

plane  (y-0)  as  assumed  by  Brown  (1978). 

2  -2 

As  the  parameter  6  ■  4kQ<h^>  increases  and  the  corresponding  value 

for  (the  wavenumber  at  which  spectral  splitting  is  assumed  to  occur) 

decreases,  it  is  necessary  to  retain  an  increasing  number  of  terms  in 

(3.18).  However,  since  the  full  wave  approach  accounts  for  specular  point 

scattering  as  well  as  Bragg  scattering  in  a  self-consistent  manner,  the 

PO 

numerical  value  for  the  total  scattering  cross  section  <a  x>  (3.15) 
should  not  depend  on  the  specific  value  of  kj,  provided  that  the  filtered 
surface  hp  satisfies  not  only  the  radii  of  curvature  criteria  but  also  the 
condition  for  deep  phase  modulation  assumed  in  reducing  (3.1)  to  the  form 
(3.15).  This  property  of  the  full  wave  solution  is  demonstrated  in  the 
next  section  where  illustrative  examples  are  presented. 

3.4  Illustrative  Examples 

In  order  to  compare  the  full  wave  solutions  for  the  scattering  cross 
sections  (3.15)  with  earlier  solutions  appearing  in  the  technical  litera¬ 
ture  (Brown  1978,  1980),  the  following  specific  form  for  the  total  surface 
height  spectrum  is  selected 


W  (v-,v~)  -  (— )S(v-,v-) 
T  x’  z  ir  x  z 


'  o  A  5  7  4 

(-)Bkq/(k  +k  )  k  <  k 

IT  —  C 

‘ 


^  0  k  >  kc  (3.23) 

where  W  is  the  spectral  notation  originally  used  by  Rice  (1953)  and  S  is 
the  notation  used  by  Brown  (1978).  For  the  assumed  isotropic  model  of  the 


ocean  surface 
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B  -  0.0046 


(3.24a) 


k2  -  v|  +  v^, (cm)-2  ,  k -  12  (cm)-1 

A  Z  C 


(3.24b) 


k  -  (335.2  V4)-Ss(cm-1),V  -  4.3  (m/S) 


(3.24c) 


in  which  V,  the  surface  wind  speed  is  given  in  m/s.  The  wavelength  for 


the  electromagnetic  wave  is 


Xq  -  2 (cm)  ,  (kQ  -  3.1416  (cm)  ) 


(3.25) 


The  mean  square  height  for  the  surface  h^  is  given  by 


2ir  kc 


JJ-t 


T(k)  kdkd$  *  f  --  -~ 
kd  kc 


(3.26) 


and  the  mean  square  slope  for  the  filtered  surface  h^  is 


2.  kd 


'/P 


k  dkd$  * 


B  r-n  ,  ,  H  +  K\ 

2  6  +  lTl  k2  J 


(3.27) 


In  (3.26)  and  (3.27)  it  is  assumed  that  k,  »  k.  The  slope  distribution 


function  is  assumed  to  be  Gaussian,  thus 


P<VV  •  -T- 

wfS 


/— 
\  4 


(3.28) 


The  directions  of  the  Incident  and  scatter  wave  normals  are  (Bahar  1981a) 


-i  .  „i  i  -  ai  -  ^  .  ni  .  Ai  - 

n  •  sin0  cos*  a  -  cos0  a  +  sin0  sind  a 
o  x  o  y  o  z 


n^  *  sin0^  coa^  a  +  cos0^  a  +  sin0^  sin<^  a 
o  y  o  y  o  z 


Thus  for  backscatter  (n*  -  -n*)0*  *  0*  *  0  •  ■  0,$^  ■  tt. 

0  0  0 


n^  “  -n1  -  -  Bin0  a  +  cos0  a 
ox  o  y 


(3.29a) 
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2k  nf 
o 


(3.29b) 


For  perfectly  conducting  surfaces  the  relative  permittivity  and  permea¬ 


bility  of  the  medium  y  <  h(x,z)  is  |er(  -*  ®,  pr  -*■  1,  thus  K^-*-  1  and 

PQ 

•>  -1  and  the  backscatter  cross  section  <a  >g(3.17)  is  given  by 


(Bahar  1981b) 


4a 
sec  0 

<0- S  ■  SPQ  —2s- 
o 

FS 


2 

tan  6 


(3.30) 


where  6nn  is  the  Kronecker  delta.  The  values  for  W  (v-,v-)  (3.22),  are 
PQ  m  x*  x 


evaluated  numerically  for  m  -  2,3.  (see  Fig.  3.3a  and  3.3b).  In  these 


plots  6  *  1.5,  and  -  VT  for  k  >  kj.  As  kj  decreases  the  number  of 
significant  terms  in  the  expression  for  <oP^>  (3.18)  increases.  The 

individual  terms  in  <0P<^>R  (i.e.,  <0P^>Rffl)  can  each  be  integrated  numer¬ 
ically  and  summed  to  obtain  <aP<^>R  or  the  integrands  of  the  individual 


terms  may  be  summed  and  numerically  integrated  once  to  give  <oP^>  . 

K 


While  the  later  procedure  is  more  efficient,  for  the  purpose  of  the 


illustrative  examples  presented  here  each  individual  term  <aPQ>  (3.19) 

Rm 


is  evaluated  separately. 


To  provide  a  basis  for  comparing  the  full  wave  solutions  with  earlier 


results,  the  normalized  backscatter  cross  sections  are  also  evaluated 


using  the  analytical  results  recently  derived  by  Brown  (1978,  1980).  Since 


his  work  is  based  on  Burrows'  perturbation  theory  (1967)  and  physical 


optics  (Beckmann  1968),  he  specifies  k.  ■  2tt/X,  (the  wavenumber  where 


spectral  splitting  is  assumed  to  occur)  on  the  basis  of  the  characteris¬ 


tics  of  the  small  scale  structure.  Thus  Brown  concludes  that  his  work 


"...clearly  demonstrates  the  merit  of  choosing  6  -  4k  <h£>  as  large  as 

possible.  On  the  other  hand,  8  must  be  less  than  one  in  order  to  satisfy 

the  basic  criteria  for  the  suitability  of  the  perturbation  technique, 

i.e.,  4kQ<h^>  <<  1".  Furthermore,  on  examining  his  illustrative  examples, 

he  notes  that  "It  would  appear  from  these  numerical  results  that  a  proper 

choice  of  kj  should  be  based  on  the  criterion  4kQ<h^>  =  0.1"...  "Any  attempt 

to  draw  a  more  positive  conclusion  about  the  choice  of  A^  would  have  to 

address  the  basic  question  of  the  dividing  lines  between  two  types  of 

scattering  mechanisms,"  Brown  goes  on  to  say,  "the  results  in  this  paper 

are  based  upon  the  assumption  that  the  scattering  is  either  physical  optics 

(or  geometric  optics  for  kQ  large  enough)  or  small  scale  diffraction. 

2  -2 

Given  this  assumption,  the  criterion  4kQ<h^>  “  0.1  seems  reasonable."  on 

the  other  hand.  Brown  (1978)  notes  that  "...the  concept  of  a  truncated  or 

filtered  spectrum  was  first  hypothesized  by  Hagfors  (1966)  in  an  attempt 

to  explain  lunar  scattering  data  and  the  observed  frequency  dependence  of 

near  normal  incidence  scattering.  More  recently  Tyler  (1976)  has 

attempted  to  definitive  Hagfors*  filter  theory  by  basing  the  spectral 

truncation  wavenumber  on  a  criterion  related  to  the  radius  of  curvature 

of  the  large  scale  surface.  Both  these  approaches  base  the  point  of 

spectral  truncation  upon  a  characteristic  of  the  large  scale  structure, 

whereas,  according  to  Brown,  "it  should  be  based  upon  the  small  scale 
2  “2 

structure  i.e.,  4kQ  "V  e<  1."  Finally,  Brown  points  out  that  his  results 
"...Indicate  a  smooth  transition  between  the  geometric  optics  and  Bragg 
scattering  regimes  which  previously  have  been  obtained  in  an  ad  hoc 


fashion." 
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Since  the  full  wave  solution,  (3.1),  accounts  for  both  physical 

optics-specular  scattering  as  well  as  Bragg  scattering  in  a  self- 

consistent  manner,  it  is  used  here  to  explore  the  basic  question  of  an 

apparent  dividing  line  between  the  two  types  of  scattering  mechanisms. 

In  Figs.  3.4a,  3.4b,  and  3.4c,  the  normalized  backscatter  cross 
PQ 

section  <a  based  on  Brown's  results  are  presented  for  B  *  0.1  (cor- 

2  -2  -1 
responding  to  oFg  “  (?£t)  *  0.0224  and  kj  -  0.95  (cm)  ).  These  figures 

show  <oPQ>-  o°Q  (the  total),  <aP(^>F  -  ^pq^o  and  <°P^>R  “  [°pq^l*  Figs. 

VV  HH 

3.4a  and  3.4b  for  <o  >  and  <cr  >  differ  slightly  from  the  results 

2 

presented  by  Brown  (1978)  for  0  0q  _<  70°  since  he  approximates  Tpp, 

P0  2 

(hpxjhp^)  (related  to  S  x)  by  its  zero  slope  approximation  Tpp,  (0,0) 

and  replaces  the  shadow  function  R  (B2  in  the  notation  used  here)  by  unity. 
These  approximations  enable  Brown  to  analytically  convert  the  two  dimen¬ 
sional  integrals  into  one  dimensional  integrals  which  he  evaluates  numeri¬ 
cally.  The  main  differences  are  near  grazing  angles  (6q  >  80°  not  shown 

in  Brown's  work)  where  shadowing  becomes  significant.  Moreover,  since 
2 

Fvh(O.O)  ■  0  he  does  not  provide  numerical  data  for  the  cross-polarized 

VH 

backscatter  cross  section  <o  >  (Fig.  3.4c).  In  view  of  Brown's  comments 

regarding  the  optimal  choice  of  6  (0.1)  and  the  sensitivity  to  his  numerical 

results  to  changes  in  8,  no  other  results  based  on  Brown's  work  are 

presented  here.  In  Figs.  3.5a,  3.5b,  and  3.5c,  the  corresponding  results 

(8  “  0.1)  based  on  the  full  wave  solutions  (3.15)  are  presented.  In  each 

of  these  figures  <cjP^>  (the  total  backscatter  cross  section)  <oPQ>c.  (3.16), 

F 

(the  cross  section  associated  with  the  filtered  surfaces)  as  veil  as 
PQ  PQ 

<0  >ri  and  <o  x>£2  (3.19)  (the  cross  sections  associated  with  h  )  are 


vv . 


Baikscatter  cross  sections  for  rough  surf ace (3. 14)  using  full  wave  formulation  (3. 15)  for  6  *  0.25 
(  -  <oPQ>  (total),  X  <aP^>„,cl  <0P^>O1 ,  4<uP^09  <aVV>. 


Fig.3.6b»  Backscatter  cross  sections  for  rough  surface(3. 14)  using  fui 
(  - <aPQ>  (total),  X«rPQ>_,  o<apQ>  A<aPQ>„,,.  <am>. 


Backscatter  cross  sections  for  rough  surface  (3.14)using  full  wave  formulation  (3.15)  for 


presented.  For  $  ■  0.1,  <0  y>gm  for  m  _>  3  are  not  shown  since  they  are  too 

P0 

small  to  affect  the  total  cross  section  <0  x>.  Furthermore,  since 

R  2  -B  22  po 

]xK(v,ns)|  -  e  -  1,  exp(-v^  <h^>)  ~  1.  and  <°r^>R2  is  ““all  compared 

PQ 

to  <o  these  results  are  very  similar  to  those  based  on  Broun  s  work. 

In  Figs.  3.6a,  3.6b,  and  3.6c  through  3.9a,  3.9b  and  3.9c,  the  scattering 

cross  sections  based  on  the  full  wave  solutions  are  presented  for  3  “  0.25, 

0.5,  1.0  and  1.5.  in  Figs.  3.9a,  3.9b,  and  3.9c,  the  term  <oP<^>Rg  is  also 

W 

presented  since  it  becomes  significant  for  3  _>  1.5.  In  Fig.  3.10,<o  >  is 

W 

presented  for  3  ■  2.0  to  demonstrate  the  importance  of  the  term  < 0  as 

3  increases  above  the  value  of  3s  1.5. 

From  the  results  based  on  the  full  wave  solutions  (Figs.  3.5  through 

3.10),  one  notices  that  while  the  individual  contributions  to  the  total 

P0 

backscatter  cross  sections  <o  are  very  sensitive  to  the  value  of  3 

(which  determines  kj)  the  value  for  the  total  cross  section  varies  slightly 

as  3  increases  from  0.1  to  1.0  and  insignificantly  for  3  ^  1.0.  This 

latter  result  is  shown  explicitly  in  Figs.  3.11a,  3.11b,  and  3.11c,  where 

only  the  total  backscatter  cross  sections  <0  >  are  plotted  for  3»0.1, 

0.25,  0.5,  0.75,  1.0,  1.25,  1.5,  1.75,  and  2.0.  It  was  not  considered 

PQ 

necessary  to  evaluate  <0  >  for  3  <  0.1  since  these  results  would  practic¬ 

ally  duplicate  those  provided  by  Brown  (1978).  In  Table  3.1,  the  com- 

2 

puted  values  for  k^  (3.26)  and  0p,g  (3.27)  corresponding  to  values  of  3  from 
0.1  to  2.0  are  listed. 

Clearly  the  above  numerical  results  indicate  that  from  the  point  of 
view  of  both  accuracy  (independence  of  the  results  on  the  choice  of  g  or 
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sections  for  rough  surface  (3.14)using  full  wave  formulation  (3.15)  for 


a  *>  (total) 
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surface  (3.l4)using  full  wave  formulation  (3.15)  for 
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Backscatter  cross  sections  for  rough  surface  (3. 14)  using  full  wave  formulation 


Backscatter  cross  sections  <a  >  for  rough  surface  (3,14)using  full  wave  formulation 


Fig.  3.11b.  Backscatter  cross  sections  for  <ory>  (total)  for  rough  surface  (3. 14) using  full  wave  formulation 

HH 

(3.15)  for  6  -  0.1,  0.25,  0.5,  0.75,  1.0,  1.25,  1.5,  1.75,  and  B  -  2.0.  <o  >. 
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kj)  and  computational  ease  the  optimal  choice  for  g  Is  1.0.  This  cor¬ 
responds  to  the  choice  of  k,  -  0.301  (cm-*)  and  X,  -  20.87  (cm)  -  10.44A  . 

a  a  o 

It  should  be  pointed  out  that  for  8  ■  1.0  not  only  does  <°^>^2  begin  to 

become  significant  (particularly  in  the  region  between  0q  *  15°  and  0q  = 

35°,  where  specular  point  scattering  and  Bragg  scattering  effects  begin  to 

R  2  3  2  2 

blend),  but  the  factors  |xvvn  )|  -  e  and  exp(-v-<h£>)  are  also  sig- 

s  y  K 

P0 

nificantly  different  from  unity.  These  factors  appearing  in  <o  >  and 

r 

PQ  PO 

<o  >  respectively  and  the  terms  <a  x>  (m  >  2)  are  not  present  In  the 
R  Rm  — 

expressions  derived  by  Brown  (1978).  The  contribution  that  the  term 
<o  >^2  makes  to  the  total  backscatter  cross  section  is  even  more  pro¬ 
nounced  for  Pj*Q  since  the  physical  optics  approximation  for  the  cross- 

PQ 

polarized  backscatter  cross  section  <0  >  is  zero.  Moreover,  since 

<oP^>F  ■  0  for  P^Q,  use  of  the  two  scale  model  of  the  composite  surface  to 
compute  the  cross-polarized  cross  sections  needs  to  be  investigated  fur¬ 
ther  (see  Section  4) . 

The  small  though  perhaps  significant  differences  in  the  full  wave 
PO 

values  for  <o  (the  total  backscatter  cross  sections) ,  as  0  increases 

from  0.1  to  1.0  are  much  smaller  than  those  obtained  from  Brown's  results. 
PO 

The  reason  why  <o  is  sensitive  to  variations  in  8  for  g  <_  1.0  is  that 

P0 

on  deriving  the  approximate  physical  optics  contribution  <0  (3.16) 

(associated  with  the  filtered  surface  h  )  from  the  full  wave  solution  (3.1) 

F 

it  was  assumed  implicitly  in  deriving  (3.30)  that  the  radii  of  curva 

ture  associated  with  the  surface  h^  satisfied  the  Kirchhoff  approximations 
for  the  surface  field  and  also  that  the  individual  specular  points  on  the 
filtered  rough  surface  h^  met  the  conditions  for  deep  phase  modulation; 
namely,  that  distances  from  the  specular  points  on  the  surface  h?  to  the 


98 


receiver  and  transmitter  were  such  that  their  contributions  were  dis¬ 
tributed  uniformly  in  phase  between  -it  and  it.  For  8  <  1.0  the  neighbor¬ 
ing  specular  points  on  h^  begin  to  merge  and  contribute  more  strongly  to 

the  backscatter  cross  section  than  when  deep  phase  modulation  applies. 

VH  HV 

For  backscatter  <a  >  -  <o  >,and  therefore  the  numerical  data  for 

HV  W 

<o  >  is  not  presented.  However,  <o  >  is  significantly  different  from 

HH  PP 

<o  >  as  0  -*■  90.  This  polarization  dependence  of  <o  >  is  missing  in 

the  physical  optics  results. 

3.5  Concluding  Remarks 

In  this  paper,  extensive  numerical  results  for  the  full  wave  back- 

scatter  cross  sections  <o  (P,Q“V  -  vertical  or  H  -  horizontal)  are 

presented.  For  convenience, and  in  order  to  compare  the  full  wave  results 

with  earlier  solutions  to  the  scattering  cross  sections,  the  full  wave 

solutions  are  expressed  in  terms  of  a  weighted  sum  of  a  scattering  cross 
P0 

section  <o  >F>  associated  with  the  filtered  surface  h^  (consisting  of 

PQ 

the  larger  scale  spectral  components)  and  a  remainder  term  <o  asso¬ 

ciated  with  the  surface  hR  (consisting  of  the  smaller  scale  spectral 
components) . 

In  an  attempt  to  draw  more  definite  conclusions  about  the  choice  of 

kj  (the  wavenumber  where  spectral  splitting  between  the  surfaces  hp  and 

2-2 

h^  is  assumed  to  occur),  the  parameter  8  -  4kQ<hR>  (that  is  related  to  kj) 
is  varied  over  a  very  wide  range  of  values  (0.1  <  6  <  2.0).  It  was  not 
necessary  to  consider  values  of  8  <  0.1,  since  these  results  would  practi¬ 
cally  duplicate  those  already  provided  by  Brown  (1978).  For  8  >  2.0  the 
numerical  evaluations  become  tedious  and  of  no  particular  value.  It  is 
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shown  that  while  there  were  small  though  significant  variations  in  the 

PQ 

values  of  <o  >  as  8  increased  from  8  *  0.1  to  6  ■  1.0  the  results  merged 
for  8  1.0.  Thus,  from  the  point  of  view  of  computational  ease  and 


accuracy,  the  choice  8  *  1.0  corresponding  to  A,  =  10.44A  is  optimal.  It 


d  o 

should  be  noted  that  in  order  to  derive  the  physical  optics  -  (specular 
point)  contribution  to  the  scattering  cross  section  (associated  with  the 
filtered  surface  h^,)  from  the  full  wave  solution,  it  was  necessary  to 
assume  that  the  radii  of  curvature  criteria  (imposed  by  the  Kirchhoff 
approximations  for  the  surface  fields),  as  well  as  the  condition  for  deep 
phase  modulation  were  satisfied.  However,  the  parameter  8  is  not  restric¬ 
ted  by  the  characteristics  of  the  small  scale  surface.  Thus  8  could  be 
assumed  to  be  much  larger  than  0.1  in  gross  violation  of  the  standard 
perturbation  criteria  without  affecting  the  results  for  the  total  scatter¬ 


ing  cross  section.  It  is  interesting  to  note  that  for  8  *  2.0  (Fig.  3.7) 


PQ  PQ 

<0  >  is  significantly  smaller  than  <o  >  for  all  values  of  0  including 

r  K  O 


near  normal  incidence.  Nevertheless,  the  computed  results  for  <o 


pQ> 


with  8  “  2.0  (the  total  backscatter  cross  sections)  are  in  agreement  with 


PQ, 


the  corresponding  results  for  0-1.0  where  <a  >  is  the  dominant  term 

r 


near  normal  incidence.  This  is  because  the  full  wave  approach  permits 
the  blending  of  specular  point  scattering  (from  h^  consisting  of  the 
larger  scale  spectral  components  of  the  surface  height)  with  Bragg  scat¬ 
tering  (from  h^  consisting  of  the  smaller  scale  spectral  components  of  the 
surface  height) .  Thus,  one  cannot  arbitrarily  neglect  the  contribution 

to  <a  associated  with  the  surface  K  ,  even  at  near  normal  incidence. 

K 

Furthermore,  since  the  physical  optics  approximation  for  the  cross- 
polarized  backscatter  cross  section  is  zero,  it  is  necessary  to  further 


A 

£ 


100 


■r.  * 

>.  • 


investigate  the  use  of  the  two  scale  model  of  the  rough  surface 
P0 

to  compute  <o  x>,  P^Q. 

Although  numerical  results  are  presented  in  this  paper  only  for 
backscatter,  the  full  wave  solution  (3.1)  is  suitable  for  the  evaluation 
of  scattering  in  arbitrary  directions.  Furthermore,  the  medium 
y  <  h(x,z)  need  not  be  perfectly  conducting  and  the  effects  of  finite 
conductivity  can  be  considered.  The  full  wave  approach  also  accounts  for 
the  scattering  of  the  lateral  and  surface  waves  (Bahar  1980a, b)  that  are 
excited  over  non-perf ectly  conducting  surfaces.  The  full  wave  approach 
can  be  applied  to  scattering  at  low  radio  wave  frequencies  as  well  as 
at  optical  frequencies. 


4 .0  Computations  of  Rough  Surface  Cross  Sections  That  Do  Not  Involve 

Spectral  Splitting 
4.1  Background 

In  Sections  2  and  3  the  full  wave  solutions  for  the  scattering 
cross  sections  were  applied  to  two-scale  models  of  composite  rough 
surfaces.  The  main  purpose  of  the  work  reported  in  these  sections 
was  to  compare  the  full  wave  solutions  for  the  like  and  cross 
polarized  scattering  cross  sections  with  the  solutions  based  on  the 
use  of  a  perturbed-physical  optics  approach.  The  question  of  the 
specification  of  k^,  the  wavenumber  where  surface  height  spectral 
splitting  is  assumed  to  occur  was  also  investigated  in  detail.  It 
was  shown  that  the  full  wave  solution  for  the  scattering  cross 
sections  can  be  expressed  as  a  weighted  sum  of  two  cross  sections.  The 
first  was  associated  with  the  filtered  surface  consisting  of  the 
large  scale  spectral  components  and  the  second  was  associated  with  the 
surface  consisting  of  the  small  scale  spectral  components  that  ride 
on  the  large  scale  surface.  Provided  that  the  large  scale  surface 
satisfied  the  radii  of  curvature  criteria  (associated  with  the 
Kirchhoff  approximations  of  the  surface  fields)  and  the  condition 
for  deep  phase  modulation,  it  was  shown  that  the  full  wave  solutions 
were  insensitive  to  the  variations  in  the  specified  value  of  k^. 

Since  the  full  wave  approach  accounts  for  specular  point 


scattering  and  Bragg  scattering  in  a  unified  self  consistent  manner, 
it  is  not  necessary  to  spectrally  decompose  the  composite  rough 


surfaces  into  two  surfaces  with  different  roughness  scales.  Thus 
in  this  section  the  scattering  cross  sections  for  all  angles  of 
incidence  are  calculated  using  the  undecomposed  form  of  the  full 
wave  solution. 

4.2  Discussion 

On  applying  the  full  wave  approach  in  Sections  2  and  3  to 
evaluate  the  like  and  cross  polarized  scattering  cross  sections  for 
two  scale  models  of  composite  rough  surfaces,  several  assumptions 
were  made  to  facilitate  the  computations.  The  first  assumption  was 
that  the  large  and  small  scale  surfaces  were  statistically  inde¬ 
pendent  (Brown  1978) .  It  would  seem  reasonable  to  make  such  an 
assumption  if  the  two  surfaces  are  results  of  independent  processes. 
This  would  be  the  case, for  example,  if  the  small  scale  roughness 
is  due  to  erosion,  while  the  large  scale  roughness  is  due  to  geo¬ 
physical  forces  that  result  in  hills  and  valleys,  or  as  in  the  case 
of  the  sea,  where  the  capillary  waves  are  dependent  on  surface  ten¬ 
sion  while  the  large  scale  rough  surface  is  generated  by  gravity 
waves.  For  the  general  case,  however,  one  cannot  assume  statistical 
independence  of  the  large  and  small  scale  surfaces. 

The  second  simplifying  assumption  that  was  made  was  that  the 
2 

mean  square  slope  n  for  the  total  surface  was  approximately  equal 

s 

2 

to  the  mean  square  slope  for  the  filtered  large  scale  surface. 

The  third  assumption  was  that  if  the  mean  square  height  of  the 


total  rough  surface  is  large  compared  to  wavelength,  the  surface 


height  characteristic  function  for  the  total  surface  is  negligibly 
small.  The  effects  of  these  simplifying  assumptions  on  the  com¬ 
puted  results  for  the  cross  sections  are  examined. 


The  physical  optics  approximation  for  the  cross  polarized  back- 
scatter  cross  section  is  zero.  As  a  result,  the  cross  polarized 
backscatter  cross  section  for  the  filtered  surface  is  set  equal 
to  zero  when  the  two  scale  model  is  used.  However,  for  backscatter 
only, the  specular  points  on  the  rough  surface  do  not  depolarize 
the  incident  wave.  Therefore,  the  justification  for  use  of  the 
two-scale  model  to  evaluate  the  cross  polarized  scattering  cross 
sections  of  composite  rough  surfaces  is  also  examined  in  this  sec¬ 
tion  of  the  report. 

4.3  Application  of  the  Full  Wave  Solution  without  Surface 
Decomposition 

The  starting  point  for  this  analysis  is  equation  (3.1)  for 
the  like  and  cross  polarized  scattering  cross  sections  of  the 
rough  surface  y  ■  h(x,z) 


«,*>>. $■  ]  <s*Wiv 

— OO  I  “  “ 

n*a 


X(vy) 


•explivxxd  +  ivzzd]dxddzd 


(4.1) 


in  which 

?d  ■  (x-x')ay  +  (z-z’)az  -  xdax  +  zdaz  (4.2) 

is  the  radius  vector  between  two  points  on  the  reference  plane  (x,z) 


The  vector  v  is 


v  *  k  (n^-n*) 
o 


va  +va  +  v  a 
xx  y  y  z  z 


(4.3) 


where  is  the  free  space  wavenumber  for  the  electromasnetic  wave 
and  n*  and  n^  are  unit  vectors  in  the  directions  of  the  incident 
and  scattered  wave  normals  respectively.  An  exp(ixot)  time  depen¬ 
dence  is  assumed  in  this  work.  The  symbol  <> denotes  the  statistical 
average  and 


—■  <SP^exp[ivy(h-h’)]> 


P2  (iif ,  n1  ]  ns  )  p  (hx ,  hz  )  dlylhz 


•X2(vy.-vy)  =  IPQ(nf,n1)x2(vy-vy) 


(4.4) 


in  which 


n(hx,hz)  is  the  unit  vector  normal  to  the  rough  surface 

f(x,y,z)  «  y-h  (x,z)  *  0  (4.5) 

Thus 

Vf  -  n|vf]  -  V(y-h(x,z))  ■  (-h  a  +a  -h  a  )  (4.6) 

ii  '  ’  xxyzz 

in  which  the  components  of  the  gradient  of  h  (x,z) 

hx  ■  3h/&x  ,  hz  «  3h/3z  (4.7) 

are  random  variables  and  p(h  ,h  )  is  the  probability  density  func- 

X  z 

tion  for  the  slopes  h  and  h  .  The  expression  for  the  scattering 

X  z 

PQ 

cross  sections  <o  x>  (4.1)  accounts  for  shadowing  and 

P2(nf,ni]n)  -  P2(nf .n1, |ng)S(nf •n)S(-ni*n)  (4.8) 

in  which  P2(n^,n*|n)  is  the  probability  that  a  point  on  the  rough 


surface  is  both  illuminated  and  visible  given  the  value  of  the 
slopes  at  the  point  (Smith  1967,  Sancer  1968)  and  P2  (n^  |  ng)  is 

its  value  at  the  specular  points  where  the  unit  vector  n  is  given 
by 


nf-n* 

l^-n1] 


v/v 


(4.9) 


The  arguments  of  the  unit  step  functions  S(-n^*n)  and  S(nf*n) 

vanish  at  points  of  the  rough  surface  where  the  incident  and 

scattered  waves  are  tangent  to  the  rough  surface.  Thus 

S(-n**n  )*■  1  and  S(n^*n  )  ■  1. 
s  s 

The  characteristic  and  joint  characteristic  functions  for  the  sur¬ 
face  height  h  are  respectively  (Beckmann  and  Spizzichino  1963) 

x(Vj)  -  <exp(ivyh)>  (4.10) 


and 


X2(vyt-vy)  -<exp[ivy(h-h')]>  (4.11) 

It  is  assumed  in  this  work  that  the  probability  density  function 
for  the  surface  height  is  jointly  Gaussian.  Thus 

x(Vy)  -  exp(-v2<h2>)  (4.12) 

and 

X2(vy*“vy)  *  exp(-v2<h2>  -  v2<hh’>)  (4.13) 

2 

where  <h  >  is  the  mean  square  height  and  <hh’>  is  the  surface 

P0 

height  autocorrelation  function.  The  coefficients  D  depend  expli 


citly  upon  the  polarization  of  the  incident  wave  (second 


superscript  Q*V  -  vertical,  Q«H  -  horizontal)  and  the  polarization 
of  the  scattered  wave  (first  superscript;  P«V,H)  the  direction  of 
the  incident  and  scattered  wave  normals  n*  and  n^  respectively, 
the  unit  vector  n  normal  to  the  rough  surface  and  the  complex  per¬ 
meability  and  permittivity  of  the  medium  of  propagation 
respectively  (Bahar  1981a,  see  Section  1).  On  deriving  (4.4)  it 
is  assumed  that  the  rough  surface  is  Gaussian  and  stationary,  thus 

the  surface  height  h  and  slopes  (h  ,h  )  are  statistically  inde- 

X  z 

pendent  (Brown  1978,  Longue t -Higgins  1957).  It  is  also  assumed 

that  for  distances  |r^|  less  than  the  surface  height  correlation 

distance,  £  ,n(h  ,h  )  -  n' (h  ’h  ').  It  has  been  shown  that  if 

the  principal  contributions  to  the  scattered  fields  come  from 

specular  points  on  the  rough  surface  (n  «  n  ),  (4.1)  reduces  to  the 

s 

physical  optics  solution  for  the  scattering  cross  section.  If, 

however,  the  roughness  scale  of  the  surface  is  small  compared  to  the 
2  2 

wavelength  (kQ<h  ><<1)  and  the  surface  slopes  hx  and  h^  are  very 
small,  (4.1)  reduces  to  the  perturbation  solution  for  the  scatter¬ 
ing  cross  sections  (Rice  1951).  Thus, in  this  case  Bragg  scatter¬ 
ing  is  accounted  for  and  the  backscatter  cross  sections  for  graz¬ 
ing  angles  are  strongly  dependent  on  polarization.  In  Section  2 
a  two  scale  model  is  adopted  to  determine  the  corresponding  full 
wave  solution  for  the  scattering  cross  sections.  To  facilitate 
the  application  of  the  two-scale  model  it  is  assumed  that  the 


small  scale  surface  h^  and  the  large  scale  filtered  surface  h^. 


are  statistically  independent  (Valenzuela  1968;  Wright  1968;  Brown 
1978).  This  assumption  is  reasonable  if  the  surfaces  h^,  and  hg 
are  results  of  independent  processes  (Brown  1978)  as  for  example, 
when  the  small  scale  roughness  is  due  to  erosion  while  the  large 
scale  roughness  is  due  to  geophysical  forces  that  result  in  hills 
and  valleys  or  as  in  the  case  of  the  sea,  where  the  capillary 
waves  are  dependent  on  surface  tension  while  the  large  scale  sur¬ 
face  is  generated  by  gravity  waves.  In  general,  however,  it 
cannot  be  assumed  that  the  large  and  small  scale  roughness  of  the 
surface  are  statistically  independent.  In  the  general  case.  If  the 
two  scale  model  is  used  to  analyze  the  problem  it  would  be  neces¬ 
sary  to  know  the  large  and  small  scale  surface  height  joint  proba¬ 
bility  density  function  for  two  adjacent  points  on  the  rough  sur¬ 
face  to  determine  (4.11)  alone. 

Since  the  full  wave  solutions  account  for  both  Bragg  scatter 
and  specular  point  scatter  in  a  unified,  self  consistent  manner, 
in  this  section  solutions  for  (4.1)  are  developed  without  adopting 
a  two-scale  model  of  the  rough  surface. 

In  Sections  2  and  3  it  has  been  noted  that  the  physical  optics 
approximation  for  the  cross  polarized  backscatter  cross  section  is 

PQ 

zero  (<o  x>j,  *  0  for  P  ^  Q). However,  even  the  large  scale  fil¬ 
tered  surface  will  depolarize  the  backscattered  field  at  non- 
specular  points  on  the  surface.  Therefore  the  present  analysis 
should  shed  more  light  on  the  evaluation  of  the  like  and  cross 


polarized  backscatter  cross  sections  and  the  suitability  of  the 

two  scale  model  even  if  it  can  be  assumed  that  the  large  and  small 

scale  surfaces  are  statistically  Independent. 

2  2  2 

Assuming  that  kQ<h  >  <<  1  and  |x|  <<  1*  the  scattering  cross 

section  (4.1)  can  be  expressed  as  follows: 


<oPQ>  -  IF(i(n£  Ix2(vy»_vy)  -  |x(v  )|  ] 


•exp[ivxxd  +  ivzzd]dxddzd 
^  IP(^(nf .n1)  Q  (nf,n-,JL) 


(4.14) 


in  which  I  ^  is  defined  by  (4.4)  and  Q  ,  the  two  dimensional 
Fourier  transform  of  (x2  —  | X I  )  depends  on  the  surface  height 


correlation  coefficient  R 


R  =  <hh'>/<h  >  . 


(4.15) 


Using  the  notation  of  Rice  (1951),  the  surface  height  spectral 


density  function  W  (v  ,v  )  is  related  to  the  two  dimensional  Fourier 

X  z 


transform  of  the  surface  height  autocorrelation  function. 
W(vx*Vz)  "  ~2  /<hh,>  expCivxXd  +  ivzZd]dxddzd 

TT  ~ 


<hh’> 


W(v  ,v  ) 
x  z 
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exp  [-iv  x,  -  iv  z,]  dv  dv 
xd  2  d  xz 


(4.16a) 


(4.16b) 


Thus  assuming  that  the  rough  surface  is  Gaussian  and  stationary,  to 

compute  the  scattering  cross  sections  (4.14)  it  is  necessary  to 

prescribe  the  two  dimensional  slope  probability  density  function 

p(h  ,h  )  (4.4)  and  the  surface  height  autocorrelation  function  or 
X  z 


Cj 
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its  Fourier  transform  (the  surface  height  spectral  density  function) 
Since  it  is  assumed  in  this  work  that  the  surface  is  isotropic. 


<hh’>  depends  only  on  the  distance  r^  *  j  rj  between  the  two  points 


(x,h,z)  and  (x' ,h* ,z' )  on  the  rough  surface.  Thus 
f  W(v  ) 

<hh’>“  2ir  I  — 7 x-z ■  -  J  (v  r,)v  dv 

1  4  o  xz  d  xz  xz 


(4.17) 


2  ,  rw<v> 

*  *  • 2  ’  J  — — 


v  dv 
xz  xz 


(4.18) 


in  which  J  is  the  Bessel  function  of  order  zero  and 
o 

2  2,2 

V  -  V  +  V 
XZ  X  z 


Since 


j  "(0)  ■  £d  J  (v  rj)/dr,]  ■  v  /2 
o  o  xz  dy  dr.*©  xz 

a 

R"(0)  -  a2/2<h2>  -  -o2/<h2>  -  -o2/<h2> 

3  X  Z 


(4.19) 


(4.20) 


(4.21) 


where 


a2  -  2ir 
s 


W(v  )  - 

xz  3  , 

— - - v  dv 

4  xz  xz 


(4.22) 


is  the  total  mean  square  slope  while  a  and  a _  are  the  mean  square 
slopes  in  the  x  and  z  directions.  Thus  for  small  values  of  the 
correlation  coefficient  is  given  by  (Beckmann  and  Spizzichino  1963, 


Brown  1978) 


R<rd>  * 1  -  - 1  -  4rd/2<h2> 


(4.23) 


where  l  is  the  correlation  distance, 
c 


'.'« rrp^;  iji*  ’.hvh  »' 
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4.4  Illustrative  Examples 

For  the  following  illustrative  examples  the  following  special 
form  of  the  surface  height  spectral  density  function  is  chosen 
(Brown  1978) 


W(v  ,v  )  -  (— )  S  (v  ,v  ) 
x*  z  it  x  z 


l(|-)Bk4/(k2+<2)4 


k  <  k 
-  c 


k  >  kc  (4.24) 


where  W  is  the  spectral  density  function  defined  by  Rice  (1951) 
and  S  is  the  corresponding  quantity  used  by  Brown  (1978).  For  the 
above  isotropic  model  of  the  ocean  surface 


B  -  0.0046 

(4.25a) 

2  2  2  -2  -1 
k  “  vx  +  vz  (cm)  »  k£  *  12  (cm) 

(4.25b) 

tc  «  (335.2V4)”*5  (cm)  1  ,  V  -  4.3  (m/s) 

(4.25c) 

In  (4.25c)  V  is  the  surface  wind  speed.  The  wavelength  for  the 
electromagnetic  wave  is 

-  3.1416  (cm)"1)  (4.26) 

Substituting  (2.24)  into  equation  (4.18)  for  the  mean  square  height 
of  the  rough  surface  yields 

2  k4 


X  *  2  (cm) 
o 


,2  B 
<h  >  =  y 


1 


1 


.2.2 
k  +k 
c 


,,2  2.2 

(kc+K  ) 


-  /.  2 ,  2.3 
3(kc+ic  ) 


- 


(4.27) 


Bk 


,  2  ,  2j.  2.3 
6k  (k  +k  ) 
c 


Thus  if  the  spectral  cut-off  point  k£  (Brown  1978)  is  much  larger 


than  k  (as  for  the  illustrative  example  (4.25)) 
B 


.2 

<h  >  « 


6<* 


r*,i 


(4.28) 
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The  surface  height  autocorrelation  function  <hh*>  (4.17)  can  be 
expressed  in  closed  form  for  k  -*■  »  .  Thus  the  surface  height 
correlation  coefficient  R(rd)  (4.15)  is  given  by  (Miller  et  al  1972) 


>■ 


9 


*  ' 
y 


i 


}■ 

> 


R(rd)  ■  [1  +  j  (<rd)2  ](Krd)K1(Kr{J)  -  (tcr^2^*^)  (4.29) 

in  which  Ko  and  are  the  modified  Bessel  functions  of  the  second 

kind  and  of  order  zero  and  one  respectively.  Since  k  »  K  and  k  >  k 

c  co 

the  above  closed  form  expression  is  used  for  R  in  this  illustrative 
example.  The  total  mean  square  slope  of  the  rough  surface  is 
obtained  on  substituting  (4.24)  into  (4.22). 


a2  -  B 
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Z n 
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,2  2 
k  +k 
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,22 

k  +k 
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11  +  l 

12  2 


,2^2 
k  -He 
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2 
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(kN-O 

k2(6tc2+15x2k2  +  Ilk4)  ^ 
c  _ c _ c 

12 (k2  +  k2)3 
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2.  2.2  ^  6  /,  2.  2V 3 


(k>0  J 
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(4.30) 


Thus  for  k  >>  k 

r  C  ,2  .  2 

2  „  fl  „  kc  +  K 

a  -  B  I—  in - = — 


-»] 


(4.31) 


For  typical  sea  surfaces  the  relative  complex  dielectric  coefficient 
at  1.5  GHz  is  given  by  (Stogryn  1971) 


e  -  42  -  139 
r 


The  slope  probability  density  function  p(hx,hz)  is  assumed  to  be 


Gaussian,  thus 

p(h  ;,h  )  - 
x  z 


exp 


it  a 


h2  +  h2  n 
x  z 


w 

In  Fig.  4.1  the  like  polarized  backscatter  cross  section  <0  >  is 

plotted  as  a  function  of  the  angle  of  incidence  0*  using  the  expres¬ 
sion  derived  in  Section  (4.3).  These  results  are  compared  with  the 
two  scale  full  wave  results  (Sections  2  and  3)  based  on  the  choice 

of  k,  (the  wavenumber  where  spectral  splitting  occurs)  correspond- 
d 

ing  to  8  *  1  (see  Section  2  equations  (2.32a),  (2.32b)  and  (2.32c)). 

W 

Both  results  yield  the  same  general  dependence  of  <0  >  on  the 

angle  of  incidence.  The  small  difference  in  level  is  primarily  due 

2 

to  the  fact  that  in  (4.3)  the  mean  square  slope  og  of  the  total 

(unfiltered)  surface  is  used, (4.31),  while  in  (2.32)  the  mean  square 

slope  aFg  for  the  filtered  surface  hp  is  used  (3.27).  It  should  be 

noted  that  in  deriving  the  expressions  for  the  scattering  cross 

sections  based  on  the  two- scale  model  (2.32),  it  was  assumed  that 
2  2 

°FS  ~  as*  ■'^1U8  results  based  on  (4.3)  are  more  accurate. 

Furthermore,  on  deriving  (2.32),  using  the  two-scale  model,  the 

quantity  x(v^)  (4.12)  is  assumed  to  be  negligible  compared  to 

X,(v  ,-v  )  (4.13)  for  r,<A  .  Since  4k  <h^>  »  3468  for  this  illus- 

trative  example,  the  resulting  approximation  is  very  good  except 

very  near  grazing  angles.  In  Figure  4.2  the  corresponding  results 

are  given  for  the  horizontally  polarized  backscatter  cross  sections 
HH 

<0  >.  It  is  interesting  to  note  that  the  full  wave  solution  (4.3) 

yields  the  proper  polarization  dependence  of  the  scattering  cross 
sections  for  all  angles  of  incidence  without  use  of  a  two  scale 
model  since  it  accounts  for  specular  point  and  Bragg  scattering  in  a 


r>.  w: 


;s  sections  <oVV>for  rough  surfaces  characterized  by  W(Vx,Vz)  g. 
model  (2.32)  and  (A)  unified  full  wave  solution  (4.4).  Relati- 


unified, self  consistent  manner.  In  Figure  4.3  the  cross  polarized 

VH  HV 

backscatter  cross  sections  <a  >  =  <o  >  are  plotted  as  functions 
of  the  angle  of  incidence.  Here  too, both  the  solutions  based  on 
the  two  scale  model  as  well  as  the  solution  derived  in  this  sec¬ 
tion  are  presented.  Unlike  the  solutions  for  the  like  polarized 

PP 

backscatter  cross  sections  <o  >  (P*V,H),  the  solutions  for  the 
cross  polarized  backscatter  cross  sections  differ  significantly, 
especially  near  normal  incidence  where  the  difference  in  level  is 
about  15db.  This  very  significant  difference  is  due  to  the  fact 
that  the  physical  optics  Approximations  for  the  cross  polarized 
backscatter  cross  section  is  zero  (Brown  1978,  see  Section  2). 

For  backscatter  the  surface  at  the  specular  points  is  normal  to  the 
incident  wave.  At  these  stationary  phase  points  no  depolarization 
occurs.  However,  since  depolarization  occurs  at  the  non  specular 
points  of  the  filtered  surface,  the  physical  optics  approximations 
for  the  cross  polarized  backscatter  cross  section  is  not  valid. 

It  is  interesting  to  note  that  for  the  two  scale  model  at  normal 
incidence 

<cjPP>/«jP^>  -  47db  (Pl*Q) 

However,  using  the  full  wave  solution  (4.3) 

<oPP>/<oPQ>  =  31db  (Pi<Q) 

The  latter  results  are  significantly  more  in  line  with  published 
experimental  results*(Long  1975). 

*See  also  NRL  report  on  "Airborne  Radar  Backscatter  Study  at  four 
frequencies",  NRL  Prob  ROZ-37,  SER:8560,  August  1966,  by  J.C. Daley 


4.5  Concluding  Remarks 

It  is  shown  in  this  section  that  two-scale  models  of  rough  surfaces 

can  be  adopted  to  obtain  solutions  for  the  like  polarized  backscatter 

cross  sections  that  are  in  reasonably  good  agreement  with  the  full  wave 

solutions  derived  in  this  section.  However,  the  two-scale  model  cannot 

be  used  to  evaluate  the  cross  polarized  backscatter  cross  sections.  The 

significant  differences  between  the  solutions  derived  in  this  section 

and  those  based  on  the  two  scale  models  (Section  2  and  3)  are  primarily 

due  to  the  fact  that  the  physical  optics  approximation  for  the  cross 

polarized  backscatter  cross  section  (associated  with  the  large  scale 

filtered  surface)  is  zero.  For  backscatter,  the  specular  points  lie 

on  portions  of  the  rough  surfaces  that  are  perpendicular  to  the  incident 

wave  normal  n^  (n  -  n^  ■  -a*-).  At  these  specular  points,  the  back- 
s 

scattered  waves  are  not  depolarized.  However  at  non-specular  points  oi 
the  rough  surface,  the  backscattered  waves  are  depolarized  (Bahar  1981b) 
Thus  it  is  important  to  note  that  even  if  a  surface  satisfies  the  radii 
of  curvature  criteria  (associated  with  the  Kirchhoff  approximations 
for  the  surface  fields) ,  the  physical  optics  approximations  for  the 
scattered  fields  may  not  be  valid  unless  for  the  given  incident  and 
scatter  angles  specular  points  exist  on  the  surface  and 
significant  contributions  to  the  scattered  fields  come  from  these 
stationary  phase-specular-points  of  the  surface.  This  explains  why  the 
physical  optics  approximations  for  the  like  polarized  backscattered 
cross  sections  are  not  suitable  for  grazing  angles  even  if  the  surface 
meets  the  radii  of  curvature  criteria  associated  with  the  Kirchhoff 


approximations . 


There  are  additional  important  reasons  for  preferring  to  use  the 
analysis  developed  in  this  section  over  those  that  are  based  on  two- 
scale  models  of  rough  surfaces.  Firstly,  if  the  two -scale  model  is 
used,  it  is  necessary  to  assume  that  the  large  and  small  scale  surfaces 
are  statistically  independent  (Brown  1978).  Secondly,  even  if  the 
assumption  of  statistical  independence  is  acceptable,  when  the  two-scale 
model  is  used,  it  is  still  necessary  to  judiciously  specify  k^  (where 
spectral  splitting  is  assumed  to  occur) .  These  problems  do  not  arise 
when  the  unified  full  wave  formulation  is  used  to  evaluate  the  scatter¬ 
ing  cross  sections. 

While  the  preliminary  results  reported  in  this  section  are  very 
significant  in  that  they  explain  why  the  earlier  solutions  based  on 
two -scale  models  of  composite  rough  surfaces  cannot  be  used  to  derive 
the  cross  polarized  backscatter  cross  sections,  more  work  needs  to  be 
done  to  apply  the  full  wave  solutions  to  more  general  (non-Gaussian, 
anisotropic) models  of  rough  surfaces. 


5.0  Concluding  Remarks  -  Summary  of  Research  During  the  Next  Six 

Months  -  Proposed  Future  Investigations. 

During  the  past  two  years  the  following  principal  contributions 
have  been  made  with  regard  to  the  application  of  the  full  wave  approach 
to  rough  surface  scattering. 

(a)  The  stationary  phase  approximation  of  the  full  wave  solu¬ 
tion  is  precisely  equal  to  the  physical  optics  solution  for  the 
scattered  fields.  Thus  the  full  wave  solution  reduces  to  the  physical 
optics  solution  when  the  expression  for  the  local  normal  to  the  rough 
surface  is  replaced  by  its  value  at  the  stationary  phase  (specular) 
points  on  the  surface.  The  questions  surrounding  the  different  exist¬ 
ing  forms  of  the  physical  optics  solutions  and  the  associated  "edge 
term"  have  also  been  resolved. 

(b)  The  full  wave  solution  reduces  to  the  geometrical  optics 
approximations  when  the  surface  integrals  are  evaluated  analytically 
using  the  steepest  descent  method. 

(c)  The  small  scale  surface  roughness  and  small  slope  approx¬ 
imation  of  the  full  wave  solution  is  precisely  equal  to  the  perturba¬ 
tion  solution  for  rough  surface  scattering.  Thus  the  full  wave  solu¬ 
tion  reduces  to  the  perturbation  solution  when  the  expression  for  the 
local  normal  to  the  rough  surface  is  replaced  by  the  normal  to  the 
reference  plane. 

(d)  Apparent  discrepancies  between  the  physical  optics  and  per¬ 
turbation  theories  have  been  resolved.  Both  theories  are  sound  if 


they  are  strictly  applied  to  surfaces  that  satisfy  the  respective 
assumed  restrictions.  Physical  optics  is  applicable  if  the  principal 


contributions  to  the  scattered  fields  come  from  specular  (station¬ 
ary  phase)  points  on  the  rough  surface.  Thus  in  order  to  apply 
the  physical  optics  appraach  to  rough  surface  scattering  it  is  not 
sufficient  to  satisfy  the  radii  of  curvature  criteria  associated 
with  the  Kirchhoff  approximations  for  the  surface  fields. 

(e)  The  full  wave  approach  provides  the  proper  limiting  forms 
for  the  far  fields  scattered  at  near  grazing  angles.  It  is  shown 
that  as  the  incident  and  scattered  fields  approach  grazing  angles 
the  plane  or  spherical  wave  approximations  for  the  far  field  are  not 
suitable.  In  this  case  on  evaluating  the  integrals  for  the  far  fields 
(in  the  wavenumber  space)  it  is  necessary  to  account  for  the  fact 
that  for  grazing  incident  and  scatter  angles,  poles  are  located  in 
the  vicinity  of  the  saddle  points.  Thus  it  is  shown  that  for  grazing 
angles  the  range  dependence  of  the  fields  is  given  by  the  error  func¬ 
tion  complement  and  the  full  wave  solutions  exhibit  the  proper  transi¬ 
tion  in  the  neighborhood  of  the  shadow  region. 

(f)  The  full  wave  approach  was  applied  to  two-scale  models  of 
rough  surfaces  and  compared  with  earlier  solutions  based  on  the  composite, 
two-scale  description  of  rough  surfaces.  The  full  wave  solutions  are 
expressed  as  a  weighted  sum  of  two  cross  sections,  the  first  associated 
with  the  filtered  surface  consisting  of  the  large  scale  spectral  com¬ 
ponents  of  the  rough  surface  and  the  second  associated  with  the  surface 
consisting  of  the  small  scale  spectral  components.  Thus  both  specular 

and  Bragg  scattering  are  accounted  for.  The  weighting  function  multi¬ 
plying  the  physical  optics  cross  section  associated  with  the  filtered 
surface  accounts  for  the  degradation  of  the  cross  section  in  the 


specular  direction  due  to  the  presence  of  the  small  scale  surface 

that  rides  on  the  filtered  surface.  It  is  shown  that  if  the 

specification  of  the  wavenumber  (where  spectral  splitting  is 

assumed  to  occur)  can  be  based  solely  on  the  characteristics  of 

the  small  scale  surface  and  if  the  mean  square  slopes  of  the  rough 

surface  are  very  small,  the  full  wave  solution  is  in  agreement  with 

solutions  based  on  a  perturbed  physical  optics  approach  (Brown,  1978) 

or  are  based  solely  on  physical  considerstions  (Wright,  1968; 

Valenzuela,  1968).  If  the  mean  square  slopes  are  not  very  small, 

Valenzuela's  solution  is  not  in  agreement  with  Brown's  solution. 

It  is  shown  through  the  use  of  the  full  wave  approach,  that  the 

difference  is  due  to  the  fact  that  in  Brown's  work  the  correlation 

distance  for  the  small  scale  surface  is  measured  in  the  reference  plane 

rather  than  along  the  large  scale  surface  as  in  Valenzuela's  work. 

Furthermore,  in  Brown's  work  the  small  scale  surface  height  is 

measured  perpendicular  to  the  reference  plane  rather  than  normal 

to  the  filtered,  large  scale  surface.  This  is  contrary  to  Burrows' 

perturbation  theory  upon  which  Brown’s  solution  is  based. 

(g)  The  controversy  between  Brown  on  the  one  hand  and  Hagfors 

and  Tyler  on  the  other  regarding  the  specification  of  the  wavenumber 

k,  where  spectral  splitting  is  assumed  to  occur  is  also  resolved 
d 

through  the  use  of  the  full  wave  approach.  As  noted  above.  Brown 
states  that  the  specification  of  k^  must  be  based  solely  on  the 
characteristics  of  the  small  scale  surface  while  Hagfors  and  Tyler 
specify  k^  on  the  basis  of  the  characteristics  of  the  large  scale 


surface.  Since  the  full  wave  approach  accounts  for  specular  point 

and  Bragg  scattering  in  a  uniform  self  consistent  manner,  k^  may 

be  specified  arbitrarily.  Provided  that  the  filtered  large  scale 

surface  satisfies  both  the  radii  of  curvature  criterion  as  well  as 

the  condition  for  deep  phase  modulation,  the  scattering  cross  sections 

based  on  the  full  wave  approach  are  shown  to  be  insensitive  to  the 

specific  choice  of  k^.  It  is  shown  that  if  one  uses  the  two -scale 

model  rough  surfaces  it  is  judicious  to  specify  k^  such  that 
2  2 

B  “  4ko<hg>  “  However  since  the  physical  optics  approximation 
for  the  cross  polarized  backscattered  cross  sections  are  zero,  the 
two-scale  model  can  only  be  used  to  evaluate  the  like  polarized  back- 
scattered  cross  sections. 

(h)  A  unified  full  wave  approach  has  been  developed  to  obtain 
the  scattering  cross  sections.  Since  this  unified  approach  does 
not  adopt  a  two-scale  model  of  the  rough  surface,  it  does  not 
artificially  separate  specular  point  scattering  from  Bragg  scattering. 
Furthermore,  using  this  uniform  approach  the  question  surrounding  the 
specification  of  the  wavenumber  k^  does  not  arise  nor  is  it  necessary 
to  assume  that  the  large  and  small  scale  surfaces  are  statistically 
independent.  It  is  shown  using  this  approach  that  the  two-scale  model 
may  be  adopted  to  evaluate  the  like  polarized  cross  sections,  but  that 
it  should  not  be  used  to  evaluate  the  cross  polarized  backscatter 
cross  sections.  The  physical  optics  cross  section  associated  with 
the  filtered  surface  vanishes  for  backscatter  because  the  specular 
points  are  located  on  regions  of  the  surface  that  are  perpendicular 


to  the  incident  wave.  Thus  the  backscattered  fields  are  not 
depolarized  from  specular  points.  For  backscatter,  depolarization 
is  due  to  scattering  from  non-specular  points.  Therefore,  as  shown 
through  the  use  of  the  full  wave  approach,  the  backscattered  fields 
are  depolarized  even  by  the  large  scale  filtered  surface- 

During  the  next  six  months  (the  final  term  of  the  current  con¬ 
tract)  the  following  phases  of  our  research  will  be  developed  further: 

(a)  Extension  of  the  preliminary  investigation  using  the  unified  full 

wave  approach  (see  Section  4). 

(b)  Investigation  of  the  dependence  of  the  scattering  cross  section 
upon  the  complex  permittivity  characterizing  the  rough  surfaces. 

(c)  Research  reported  in  the  Interim  Technical  Report  on  "Scattering 
Cross  Sections  for  Composite  Models  of  Non-Gaussian  Rough  Surfaces  for 
Which  Decorrelation  Implies  Statistical  Independence"  revised  and  re¬ 
submitted  for  publication  in  scientific /technical  journals. 

(d)  Computer  programs  updated  to  reflect  recent  advances  in  the 
analytical-numerical  work.  A  detailed  listing  of  the  computer  programs 
will  be  presented  in  the  final  report.  The  documentation  will  also 
include  flow  charts  and  relevant  comments  for  the  convenience  of  the 
user. 

It  is  proposed  that  the  following  topics  be  considered  for 
future  investigation  as  an  extension  of  the  current  contract: 

(a)  Scattering  and  depolarization  due  to  rough  surfaces  covered  by 
vegetation 

(b)  Scattering  due  to  variations  in  the  complex  permittivity,  e(r), 
that  characterizes  the  rough  surface  (mixed  path  propagation  using 
the  full  wave  approach) 


(c)  Application  to  synthetic  aperture  radars 

(d)  The  dependence  of  the  scattering  cross  sections  on  different 
rough  surface  height  spectral  density  functions 

(e)  Extension  of  the  research  on  scattering  by  non-Gaussian  rough 
surfaces 

(f)  Application  of  the  full  wave  approach  to  anisotropic  rough 
surfaces. 

The  extension  of  the  current  research  could  be  conducted  over 
a  30-month  period  at  approximately  the  current  level  of  support.  The 
principal  investigator  will  consider  additional  applications  that  are 
of  interest  to  the  contractor. 
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